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An  analysis  is  developed  for  predicting  the  near  zone  fields  scattered  by 
perfectly-conducting,  open-ended,  semi -inf i nite  circular  and  rectangular  waveguides  (or 
ducts)  which  are  excited  by  an  external  electromagnetic  (EM)  plane  wave.  These  waveguides 
are  terminated  inside  by  various  structures  such  as  a  planar  impedance  surface,  or  a 
disk/blade  structure  with  and  without  a  hemispherical  or  a  conical  hub.  The  analysis  is 
based  on  a  combination  of  asymptotic  high  frequency  techniques  such  as  ray  methods,  the 
equivalent  current  approach,  and  the  physical  theory  of  diffraction,  with  the  usual  modal 
techniques,  all  of  which  are  used  in  conjunction  with  the  multiple  scattering  method 
(MSM).  Basically,  the  scattering  matrices  in  the  MSM,  which  account  for  the  multiple  wave 
interactions  between  the  open  end  and  the  interior  termination  are  found  efficiently  via 
the  asymptotic  high  frequency  techniques.  The  latter  techniques  provide  relatively  simple 
expresssions  not  only  for  the  fields  scattered  from  the  rim  edge  at  the  open-end  of  the 
ducts,  but  also  for  the  amplitudes  of  the  modes  coupled  into  the  duct  by  the  incident 
plane  wave  as  well  as  the  reflection  coefficients  of  the  modes  at  the  open  end.  Some 
numerical  results  are  provided  for  the  near  and  far  field  bistatic  scattering.-.  Whenever 
possible,  a  comparison  via  other  approaches  is  provided  to  illustrate  the  accuracy  of  the 
numerical  results  for  the  elements  of  the  scattering  matrices  of  the  MSM,  since  the 
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accuracy  of  these  elements  ultimately  decide  the  accuracy  of  the  complete 
solution  of  the  scattering  by  the  inlets  as  developed  in  this  report.  For 
the  circular  waveguide  case,  the  accuracies  of  extensive  calculated 
backscattering  returns  have  been  substantiated  by  measurements.  The  method 
is  also  applied  to  evaluating  plane  wave  scattering  by  rectangular  and 
semi-circular  waveguides  mounted  on  a  ground  plane. 
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I.  INTRODUCTION 


In  the  present  report,  an  efficient  analysis  is  developed  for 
predicting  the  near  zone  fields  scattered  by  perfectly-conducting, 
open-ended,  semi -infinite  circular  and  rectangular  waveguides  and  also 
semi -infinite  semi -ci rcular  and  rectangular  waveguides  mounted  on  an 
infinite  ground  plane,  which  are  excited  by  an  external  electromagnetic 
(EM)  plane  wave.  These  waveguides  which  model  simple  inlet 
configurations  are  terminated  inside  by  various  structures  such  as  a 
planar  impedance  surface,  a  disk,  a  disk-blade  geometry,  and  also  by 
hemispherical  and  conical  hubs  on  disks  and  disk-blade  geometries  as 
shown  in  Figures  1(a)  through  1(g). 

The  present  analysis  employs  a  combination  of  high  frequency 
techniques  such  as  the  geometrical  theory  of  diffraction  (GTD)  [1,2,3] 
and  its  modifications  based  on  the  equivalent  current  method  (ECM) 
[4,5,6],  as  well  as  the  physical  theory  of  diffraction  (PTD)  [6,7],  all 
of  which  are  used  in  conjunction  with  the  self-consistent  multiple 
scattering  method  (MSM)  [8],  as  indicated  in  a  previous  analysis 
[9,10]. 

As  is  well  known,  the  GTD  employs  rays  to  describe  the  phenomenon 
of  wave  radiation,  propagation,  scattering  and  diffraction  at  high 
frequencies.  Furthermore,  the  GTD  is,  in  general,  found  to  remain 
accurate  even  for  moderately  high  frequencies.  In  the  present 
analysis,  the  GTD  together  with  its  modification  at  caustics  and  a 
confluence  of  caustics  and  shadow  boundaries,  which  are  based  on  the 
ECM  as  well  as  the  PTD,  is  employed  to  calculate  the  elements  of  the 
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Figure  1(a).  EM  scattering  by  a  circular  Inlet  with  a  planar  impedance 
termination  Inside. 


Figure  1(b).  EM  scattering  by  a  rectangular  Inlet  with  a  planar 
Impedance  termination  Inside. 
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Figure  1(c).  EM  scattering  by  a  circular  inlet  with  a  conical  hub. 


ie\ 


Figure  1(d).  EM  scattering  by  a  circular  inlet  with  a  hemispherical 
hub. 
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FRONT  VIEW: 


HUB  WITH  HUB  WITH 

BLADE  STRUCTURE  SHORT  END 

Figure  1(e).  Different  planar  terminations  on  which  the  hub  is 
attached. 


Figure  1(f).  EM  scattering  by  a  rectangular  waveguide  mounted  on  a 
ground  plane. 

PLANAR  IMPEDANCE  ,  HUB, 


Figure  1(g).  EM  scattering  by  a  semi-circular  waveguide  mounted  on  a 
ground  plane. 


junction  scattering  matrices  which  are  required  in  the  MSM.  These 
scattering  matrices  characterize  the  scattering  properties  of  the 
junctions  (T)  and  (n)  at  z=0  and  z=-L,  respecti vely ,  which  are 
indicated  in  Figures  2(a)  and  2(b).  The  multiple  wave  interactions 
between  junctions  (T)  and  (n)  is  accounted  for  in  a  self  consistent 
fashion  in  the  MSM  through  the  use  of  the  junction  scattering  matrices. 
It  is  noted  that  the  self-consistent  MSM  procedure  used  here  is 
essentially  the  same  as  the  generalized  scattering  matrix  technique 
(GSMT)  [11]. 

The  format  of  this  report  is  as  follows.  First,  the  self- 
consistent  MSM  formulation  is  discussed  in  Section  II  where  the 
junction  scattering  matrices  are  defined.  Next,  the  GTD/ECM/PTD  based 
analysis  to  determine  the  elements  of  these  junction  scattering 
matrices  in  a  relatively  efficient  manner  is  developed  in  Section  III. 
Some  numerical  results  illustrating  the  behavior  of  the  junction 
scattering  matrices  are  also  presented  in  Section  III.  Finally,  these 
relatively  simple  expressions  for  the  elements  of  the  junction 
scattering  matrices  are  combined  via  the  MSM  procedure  to  obtain 
numerical  results  for  the  near  field  scattered  by  the  circular, 
semi-circular  and  rectangular  inlet  configurations  of  Figure  1;  several 
such  numerical  results  and  comparisons  with  measured  ones  are  presented 
in  Section  IV  together  with  a  discussion  of  these  results.  An  e+JU)t 
time  convention  for  the  electromagnetic  fields  is  assumed  and 
suppressed  throughout  this  report. 
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II.  SELF-CONSISTENT  MSM  FORMULATION 


As  mentioned  in  Section  I,  the  method  of  analysis  which  will  be 
employed  here  to  predict  the  near  field  bistatic  EM  scattering  by  the 
configurations  in  Figure  1  is  based  primarily  on  the  GTD  and  the  ECM, 
as  well  as  the  PTD,  together  with  the  self-consistent  MSM.  In  the  MSM 
based  analysis,  the  basic  scattering  mechanisms  are  isolated  and 
identified  as  being  associated  with  the  scattering  by  junctions  (T)  and 
^i)  as  in  Figures  2(a)  and  2(b).  These  junctions  can  likewise  be 
identified  in  open-ended  semi-circular  and  rectangular  waveguides  on  a 
ground  plane.  The  multiple  scattering  between  the  junctions  is 
calculated  via  a  self-consistent  procedure.  Such  a  procedure  requires 
a  knowledge  of  the  generalized  scattering  matrices  [Sji],  [S^],  [$2lL 
[S22],  and  [Sr]  for  the  junctions  (7)  and  (il)  .  These  generalized 
scattering  matrices  [9]  are  directly  associated  with  the  canonical 
scattering  events  shown  in  Figures  3,  4,  5  and  6  for  an  open-ended 
semi -infinite  waveguide  with  an  impedance  termination  at  z=-L.  The 
waveguide  cross-section  is  shown  to  be  arbitrary  In  Figures  3,  4,  5  and 
6;  however,  as  stated  previously,  only  circular  and  rectangular 
cross-sectional  waveguides  are  of  interest  in  the  present  study. 

Besides  an  impedance  termination,  other  terminations  studied  include  a 
disk,  a  disk-blade,  and  hemispherical  and  conical  hubs  on  disk  and 
disk -blade  structures.  The  concept  of  the  scattering  matrix  in 
microwave  circuit  analysis  Is  based  on  the  propagating  modes  within  the 
circuit  and  It  characterizes  waveguide  discontinuities  or  circuit 
properties  of  a  microwave  network.  Thus,  the  circuit  scattering  matrix 
for  microwave  networks  is  restricted  to  interior  (guided  wave)  regions. 
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On  the  other  hand,  the  polarization  scattering  matrix  is  defined  for 
exterior  regions  [12].  The  generalized  scattering  matrix  extends  the 
concept  of  the  microwave  circuit  scattering  matrix  to  include 
evanescent  modes  as  well  as  to  include  the  polarization  scattering 
matrix  for  exterior  scattering.  Thus,  the  generalized  scattering 
matrix  approach  is  useful  for  solving  scattering  problems  which  involve 
a  coupling  between  interior  and  exterior  regions  as  in  the  problems  of 


Figure  1.  The  elements  of  the  generalized  scattering  matrices 
associated  with  the  junctions  (T)  and  (7^  as  in  Figures  3,  4,  5,  and  6 


are  defined  next. 

Let  E1  denote  the  electric  field  intensity  of  the  incident 
electromagnetic  plane  wave.  One  may  express  E1  as: 


(1) 

(2a) 


E 


i 


.ri  - 
-jk  t 

e 


where 

ic1  =  -Msin^cos^x  +  sine^in^y  +  cos^z) 
and 


(2b) 


(3) 


A  A  A 

r  -  xx  +  yy  +  zz 


(4) 


8 


Figure  4.  Scattering  by  only  the  rim  at  the  open  end  of  the  inlet. 
Scattering  occurs  both  in  the  exterior  as  well  as  the 
interior  regions  of  the  inlet. 
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Figure  5. 


Radiation  and  reflection  at  the  open  end  due  to  a  waveguide 
mode  incident  at  the  opening. 


Figure  6.  Reflection  of  waveguide  modes  from  the  Impedance  surface 
termination  at  z=-L. 
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It  is  noted  that  k  refers  to  the  free  space  wave  number  and  r  refers  to 
the  position  vector  of  an  observation  point  at  P. 

One  may  now  describe  the  scattering  matrix  [S^]  (see  Figures  3 

—  c  r» 

and  4).  The  scattering  matrix  [S^]  relates  the  electric  field  E 
scattered  from  only  the  open  end  at  z=0  to  the  field  E1  which  is 
incident  on  the  open  end  as  follows: 


£SO 

sSO 

sSO 

X 

X0 

x<t> 

Es° 

s5° 

Ss° 

y 

yQ 

y$ 

^so 

sSO 

sSO 

z 

ze 

Z4> 

_  _ . 

_ 

— 

where 


(5) 


(P)  •  E^°x  +  Ej°y 


*  zp 


(6) 


and 


csn] 


-so 

sxe 

-so 

ye 

-so 


5S° 

X<t> 

5-SO 

y4» 

r-SO 


ze  z* 


(7) 


The  scattering  matrix  [Sj^]  converts  the  waveguide  modal  fields 
incident  at  the  open  end  (z=0)  from  the  interior  region  (z<0)  to  the 

fields  radiated  by  these  modes  from  the  open  end  as  in  Figures  3 
and  5.  The  modal  electric  field  E*  within  the  waveguide  region  (z<0) 

may  be  represented  in  the  usual  manner  [13]: 
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E+  =  H  [A+  e  e  '  nm  +  B*  (e  ±  ze  )  e  ^  nm  ] 
g  ““  L  nm  nm  nm  nm  znm  J 


(8) 


Here,  e  and  e  denote  the  transverse  (to  z)  electric  vector  mode 
nm  nm 

functions  for  the  TE  and  TM  modes,  respectively  (with  TE  and  TM 

nm  nm  r  J  ' 

I 

being  defined  with  respect  to  z).  Likewise,  the  3  and  3  denote 

nm  nm 

the  propagation  constants  of  the  TE  and  TM  modes,  respectively. 

nm  nm 

Also,  n  and  m  refer  to  the  modal  indices  associated  with  the  transverse 
eigenvalues  (in  x  and  y  for  the  rectangular  waveguide  case,  or  in  <j>  and 
p  for  the  circular  waveguide  case)  associated  with  the  waveguide 

A 

eigenfunctions  or  modes.  It  is  noted  that  e  Is  the  z-component  of 

znm 

I 

the  TM  modal  electric  field;  clearly,  e  for  the  corresponding 
nm  znm  r  3 

TEnm  case  vanishes  by  definition.  The  superscripts  +  and  -  in  (8) 


refer  to  modes  propagating  in  the  +z  and  -z  directions,  respectively. 
It  is  convenient  to  define  the  magnetic  field  H*  in  the  waveguide 
region  following  the  representation  for  the  electric  field  E~  in  (8); 
thus  [13], 


-+  *  Z  +  jBZ 

H  =  II  [A4  C±h  +  h  z)e  nm  +  B*  (±R  )e  nm  ] 

g  L  nnr  nm  znm  1  nm'  nm  J  • 

3  mn  ( g ) 

The  A*  and  B*  in  (8)  and  (9)  are  the  modal  coefficients.  If  Er0(P) 
nm  nm 

denotes  the  electric  field  at  P  exterior  to  the  waveguide  region,  which 
Is  radiated  by  the  modes  that  Impinge  on  the  open  end,  then  the 
scatering  matrix  [S^  relates  Ero  to  E^  as  follows: 
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(10) 


(11) 


(12) 


(13a) 


(13b) 


The  scattering  matrix  [S,,^]  describes  the  transformation  or  the 

coupling  of  the  incident  plane  wave  field  into  the  waveguide  modes  as 

Illustrated  In  Figure  4.  It  Is  clear  that  can  therefore  relate 

E"  to  V  as 
9 
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where 


[A-  ] 

L  nmJ 

[B“  ] 

L  nmJ 

[s 

fS21^ 

^S9mJ  ^S<j>nm- 
^S9ntJ  ^niJ 


<}>nmJ 

5  1 
4>nmJ 


(14) 


(15) 


It  is  noted  that  the  problem  of  determining  [S21")  is  related  to  the 
problem  of  determining  [S12I  via  reciprocity  if  P  is  moved  to  the  far 
zone  in  the  reciprocal  (or  [S12D  problem  tc  correspond  to  a  local 
plane  wave  incidence  in  the  original  (or  [S21D  problem.  The  precise 
relationship  between  [S21 3  and  [S12]  will  be  discussed  in  Section 
1 1 1  -B . 

The  scattering  matrix  CS22 1  is  a  modal  reflection  coefficient 

matrix  which  is  associated  with  the  intei action  illustrated  in  Figure 

5.  In  particular,  the  elements  of  [S22l  describe  the  reflection 

coefficients  associated  with  the  TE  and  TM  modes  reflected  back 

nm  nm 

from  the  open  end  (at  z=0)  into  the  waveguide  region  (z<0)  when  either 
a  TE  or  a  TM  mode  is  incident  on  that  open  end  from  within  the 

pq  pq 

waveguide.  Thus,  the  matrix  CS 22 re^a^es  ^  to  ^  as  follows: 


f  R^  1 
1  nm;pqJ 

[Rhe  1 

1  nm;pqJ  | 

X,i 

__ 

[Reh  I 

1  nm;pq] 

[Re3  1  ' 

L  nm;pqi  | 

[»PV 

where 
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cs_] .  ~[CPq]  tR™1Pqi' 

foeh  i  rRee  i 
L  nm;pqJ  L  nm;pqJ 


The  meaning  of  R  for  example.  In  (17)  Is  the  following.  A  TM 

nm;pq  pq 

(or  e  type)  mode  with  modal  amplitude  Bpq  which  is  incident  at  the  open 

end  is  partly  transformed  (or  coupled)  into  a  reflected  TE  (or  h 

nm 

type)  mode  with  modal  amplitude  A"  =  Rhe  B+  . 

nm  nm;pq  pq 

The  scattering  matrix  [Sj*],  like  [S22L  is  also  a  reflection 
coefficient  type  matrix  which  is  associated  with  the  termination  at 
z=-L  (see  Figure  5).  Thus,  one  may  write: 


10 

jO_ 

~frhh  1  frhe  1 
L  nm;pqJ  >•  nm;pqJ 

[reh  ]  [ree  1 
L  nm;pqJ  Llqm;pqJ 


tApq) 

[Bpq] 


where 


tv- 


[rhh  1  frhe  1 

L  nm;pqJ  1  nm;pqJ 

[reh  ]  [ree  ] 

L  nm;pqJ  L  nm;pqJ 


As  before. 


rnm" nn  re1ates  A*  to  B“  via  A* 
nm;pq  nm  pq  nm 


nm;pq 


At  any  given  operating  frequency,  the  waveguide  region  can  support 


a  finite  number  of  propagating  modes  and  an  infinite  number  of 
evanescent  (non-propagating)  modes.  Therefore,  the  matrices  [S^l* 
[S21 3 •  CS22^  and  [Sr]  are  of  infinite  order  to  Include  the  Infinite 
number  of  evanescent  modes.  However,  even  though  the  matrices  [S12L 
[S21L  [S22]  and  [Sr]  are  of  infinite  order  in  a  formal  sense,  one 
needs  to  retain  only  a  finite  number  of  the  elements  of  these 
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scattering  matrices  in  practice,  because  the  distance  "L"  shown  in 
Figure  1  is  generally  large  enough  to  where  the  infinite  number  of 
evanescent  modes  generated  at  junction  ^T)  contribute  negligibly  at 
junction  (n)  ,  and  vice  versa.  The  finite  number  of  elements  of  the 
scattering  matrices  which  are  retained  in  practice  thus  correspond  to 
only  the  finite  number  of  all  the  propagating  (or  non  evanescent)  modes 
which  can  exist  within  the  waveguide  region.  If  the  distance  L  in 
Figure  1  is  small  enough  so  that  the  "lower  order"  evanescent  modes 
become  important,  then  one  must  include  these  modes  but  one  can  still 
ignore  the  contribution  from  all  the  "higher  order"  evanescent  modes 
when  the  latter  is  negligible.  In  either  case,  one  always  retains  a 
finite  number  of  elements  in  the  scattering  matrices  which  may 
otherwise  formally  be  of  infinite  order.  The  number  of  evanescent 
modes  which  need  to  be  included  in  the  special  case  when  L  is  very 
small  is  dictated  by  how  many  of  these  are  needed  for  convergence;  it 
has  been  known  from  experience  that  even  when  L  *  0,  the  inclusion  of 
only  three  to  five  evanescent  modes  are  sufficient  for  convergence. 

It  is  observed  from  Figure  1  that  one  may  describe  the  field 
scattered  by  the  inlet  as  consisting  of  two  parts.  The  entire 
scattered  field  thus  consists  of  the  contribution  from  just  the  open 
end  by  itself,  together  with  the  field  radiated  from  the  open  end  as  a 
result  of  all  the  interactions  between  the  open  end  and  the  termination 
that  arises  from  the  incident  field  which  is  coupled  into  the  waveguide 
region.  Referring  to  Figure  7,  it  is  observed  that  the  incident 
electric  field  [E1]  =  [a9]  at  the  open  end  is  scattered  by  junction  Cl) 
to  produce  a  scattered  field  [Esol  given  by  (5)  as 
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[ES°]  -  [Sn][Ei] 


(20) 


where 


[  Eso  ]  = 


:so 

'x 

-so 


-so 


(21) 


A  part  of  the  Incident  field  is  coupled  at  ^7)  into  the  waveguide 
region;  this  field  then  becomes  incident  at  junction  ^7)  from  which  it 
is  subsequently  reflected.  The  field  reflected  from  junction  is 

incident  back  at  where  it  undergoes  further  scattering  into  the 
exterior  and  interior  regions,  and  so  on,  thus  giving  rise  to  multiple 
wave  interactions  between  junctions  and  .  The  fields 
resulting  from  these  multiple  interactions  may  be  expressed  in  a 
convergent  Neumann  series  as  done  by  Pace  and  Mittra  [14];  however,  an 
alternate  procedure  based  on  a  self-consistent  method  leads  directly  to 
the  same  result.  The  latter  self-consistent  MSM  based  analysis  will  be 
employed  here.  Let  [E12]  represent  the  total  field  Incident  at  (T) 
from  (n)  after  taking  all  the  multiple  Interactions  into  account. 
Likewise,  let  [E21]  represent  the  total  field  incident  at  (n)  from 
(T)  after  taking  all  the  multiple  interactions  Into  consideration. 
Then,  the  total  scattered  field  in  the  exterior  region  denoted  by  [Es] 
consists  of  a  superposition  of  [Eso]  and  [Em0],  where  [Emo]  is  the 
field  scattered  Into  the  exterior  region  when  [E^]  Is  incident  on  (7). 
Therefore, 
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[ES]  =  [ES0]  +  [E010] 


(22) 


where  [Em0]  may  be  expressed  as 
[Em°]  =  [ S 12 ]  [E12] 

The  expressions  for  [E^  anc*  ^21^  may  be  written  in  terms  of  CS21  ] 
and  [S22];  namely, 

[E21]  -  [S2  IfE1  ]  *  [S,J[E .,]  (24) 

z=0  z=0 

and 

[E,?]  ■  [P]  [S  ][P][E21]  (25) 

z=0  z=0 

in  which  [Sr]  denotes  the  reflection  coefficient  scattering  matrix 
pertaining  to  junction  (n)  as  defined  in  (19),  and  [P]  is  a  diagonal 
matrix  accounting  for  the  phase  delay  along  the  length  L  between 
junctions  (T)  and  (fj)  .  Eliminating  [E]^  between  (24)  and  (25) 
yields 

([I]  -  [S22][P][S.][P])  [E  ]  .  [s21] [E1 ]  ,  (26) 

z=0 

where  [I]  is  an  infinite  order  identity  matrix.  Hence,  it  follows 
that 

[E21]  -  ([l]-[S22](P](Sr]{P])-1[S21][Ei]  .  (27) 

za0 
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Figure  7. 


Multiple  Interactions  between  junctions  CD  and  (i^  for 
calculating  the  total  backscattered  field. 
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From  (25)  and  (27),  [E12I 


z=0 


becomes 


[P]tSr][P]([l]-[S22][P][Sr][P])-1[S21][E1] 


(28) 


Incorporating  the  above  result  into  (23)  yields 

[Em0]  =  [S12][P][Sr3[P]([n  -  CS22][P][Sr][P])’1tS21][E1]  .  (29) 


Finally,  combining  (20)  and  (29)  according  to  (22)  leads  to  the 
self-consistent  expression  for  the  total  scattered  field  [Es]  as 

[ES]  =  {[Su]+[S12][P][Sr][P]([I]-[S22][P][Sr][P])-1[S21]}[Ei]  . 

(30) 

The  scattering  matrices  [$n],  [S21L  and  [S22]  can  be  found  in 

a  relatively  simple  manner  via  asymptotic  high  frequency  techniques 
(such  as  GTD,  UTD,  ECM  or  PTD),  as  shown  in  the  next  section. 

Before  proceeding  further  with  the  development  of  the  expressions 
for  the  elements  of  the  scattering  matrices,  it  is  important  to  mention 
that  these  scattering  matrices  are  somewhat  different  from  the 
conventional  scattering  matrices  due  to  the  fact  that  they  contain 
range  information  (because  they  are  valid  for  field  calculations  in  the 
near  zone  of  the  Inlet  opening).  Secondly,  one  notes  that  the  present 
analysis  is  restricted  to  angles  of  incidence  and  scattering  which  are 
about  60°  or  less  from  the  forward  axis  of  the  Inlets. 
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In  the  far  zone  calculations,  it  is  more  common  to  define  the 
scattering  matrices  so  that  they  are  independent  of  the  range  (i.e., 
the  far  zone  distance  from  the  origin  of  the  coordinate  system  shown  in 
Figure  1).  The  alternative,  more  common  definitions  for  the 
generalized  scattering  matrices  used  in  far  zone  calculations  are 
included  in  Appendix  V  for  the  sake  of  completeness  and  convenience. 

It  is  noted  that  the  results  calculated  from  (30)  in  the  limit  as  the 
near  zone  point  is  allowed  to  recede  to  the  far  zone  should  reduce  to 
those  obtained  from  the  far  zone  scattering  described  in  Appendix  V. 

III.  ANALYSIS  FOR  THE  ELEMENTS  OF  THE  SCATTERING  MATRICES 

A.  Development  of  [S|j] 

The  scattering  matrix  [Sn]  describes  the  field  scattered  from 
only  the  open  end  (at  z=0)  of  the  geometries  in  Figure  1  when  these 
geometries  are  illuminated  by  an  external  source.  In  the  present 
study,  the  illumination  is  an  external  EM  plane  wave  and  the 
observation  point  P  is  also  located  in  the  exterior  region; 
furthermore,  P  could  lie  either  In  the  near  zone  or  the  far  zone  of  the 
aperture  (at  z=0)  corresponding  to  the  open  end  of  the  inlet 
configurations  in  Figure  1.  According  to  the  ray  picture  afforded  by 
the  GTD,  the  fields  scattered  by  just  the  open  end  arise  only  from  the 
diffraction  of  the  incident  wave  by  the  rim  corners  and  edges  (at  the 
inlet  opening).  Thus,  let  Eso  and  Hso  denote  the  electric  and  magnetic 
fields,  respectively,  which  are  singly  diffracted  by  the  inlets  in 
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Figure  1;  in  addition,  multiply  edge  diffracted  fields  (Esm,  Hsm), 
which  are  produced  by  rays  that  undergo  multiple  diffractions  across 
the  aperture,  also  contribute  to  the  scattered  field  at  P.  However,  in 
the  present  study,  only  Eso  and  Hso  are  included;  whereas,  Esm  and  Hsm 
are  not  considered.  These  multiple  interactions  may  become  important 
if  the  aperture  dimensions  are  not  sufficiently  large  in  terms  of  the 
wavelength. 

Actually,  Eso  (and  Hso)  is  the  sum  of  the  fields  of  all  the  singly 
diffracted  rays  which  are  initiated  at  various  points  on  the  rim  edge 
by  the  incident  field  in  accordance  with  the  law  of  edge  diffraction, 
which  is  a  consequence  of  one  of  the  basic  postulates  of  the  GTD.  That 
postulate  deals  with  Keller's  generalization  of  Fermat's  principle  to 
include  diffracted  rays.  In  the  case  of  the  rectangular  inlet  opening, 
the  rays  diffracted  from  the  corners  of  the  inlet  rim  also  contribute 
to  the  diffracted  fields  Eso  and  Hso  in  addition  to  the  edge  diffracted 
rays.  The  points  of  edge  diffraction  can,  in  general,  migrate  around 
the  rim  as  a  function  of  the  observation  point  and  incident  field 
direction.  In  the  case  of  the  circular  inlet  the  entire  rim  edge  or 
portions  of  the  rim  can  contribute  to  the  edge  diffracted  field  for 
certain  observation  points;  in  such  instances  there  exists  a  caustic  of 
the  edge  diffracted  rays  at  or  near  P  and  the  GTD  fails  because  of  the 
continuum  of  rays  that  must  now  reach  P  in  this  case.  Also,  away  from 
the  caustics  of  diffracted  rays  there  can,  in  general,  be  four  points 
of  diffraction  on  the  circular  rim  which  contribute  to  the  field  Eso 
when  it  is  evaluated  for  a  certain  set  of  oberservation  points  P; 
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whereas,  at  other  observation  points  P  only  two  points  of  edge 
diffraction  on  the  circular  rim  contribute  to  Eso.  Again,  a  direct 
application  of  the  GTD  fails  in  the  range  of  observation  points  within 
which  such  a  disappearance  of  four  points  into  two  points  of 
diffraction  occurs  on  the  rim  edge.  Typically,  such  a  disappearance  of 
rays  results  when  three  of  the  four  points  merge  together.  Conversely, 
one  can  also  experience  a  transition  from  two  to  four  rays  diffracted 
to  P.  In  order  to  overcome  some  of  the  above  problems  in  the  GTD  which 
occur  for  a  certain  range  of  aspects,  it  is  necessary  to  use  the 
equivalent  current  method  (ECM).  In  the  case  of  the  rectangular  inlet, 
the  use  of  ECf<  also  automatically  takes  into  account  the  presence  of 
the  corners  at  the  inlet  opening.  It  is  noted,  of  course,  that  the 
effect  of  the  corners  is  taken  into  consideration  in  an  approximate 
fashion  in  the  ECM;  such  an  approximation  is  also  present  if  one 
employs  the  GTD  approach  alone  since  the  GTD  corner  diffraction 
coefficient  that  is  available  at  this  time  is  based  partly  on  heuristic 
arguments  [3]  and  is  thus  approximate.  In  regions  where  GTD  is  valid, 
the  ECM  generally  blends  into  the  GTD  solution  provided  the  inlet 
opening  is  sufficiently  large  in  terms  of  the  wavelength.  Thus,  one 
could  switch  over  from  the  ECM  based  solution  to  a  GTD  solution  away 
from  regions  where  GTD  fails.  Such  a  switch  would  yield  a  numerically 
more  efficient  approach  for  evaluating  the  elements  of  [Sn]  because 
the  ECM  involves  an  integration  around  the  rim  edge  (of  the  inlet), 
whereas  the  GTD  involves  no  integration.  In  the  present  work,  this 
kind  of  switching  has  not  been  used  for  the  calculation  of  the  [S^] 
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matrix.  The  ECM  method  is  employed  at  all  aspects  of  interest  in  this 
study. 

In  the  ECM  [3-6],  equivalent  currents  Ieq  and  Meq  of  the  electric 
and  magnetic  type,  respectively,  are  located  at  the  rim  edge  as  shown 
in  Figure  8  for  an  arbitrarily  shaped  inlet  opening.  These  equivalent 
rim  currents  radiate  in  "free  space"  to  produce  the  scattered  field. 

As  indicated  in  [3,6],  Ieq  and  Meq  are  actually  deduced  from  GTD. 

Thus,  even  though  GTD  itself  is  not  valid  within  diffracted  ray  caustic 
regions,  it  can  still  indirectly  provide  the  strengths  of  the 
equivalent  currents  which  radiate  fields  that  remain  valid  within  those 
caustic  regions  in  the  ECM  based  approach.  From  [6], 


I 

eq 


(*') 


DsOM^go.g) 
|/sinSQsine  |  ^  'Z  ^ 


and 


Dh(^»,;g0.g)  - 

|/sin6Qsin6  |  ^ 


(31) 


(32) 


Z0  =  free  space  impedance,  and 
Y0  =  1/Z0  . 

One  can  incorporate  Ieq  and  Meq  of  (31)  and  (32)  into  the  radiation 
integral  for  calculating  the  scattered  field  Eso;  namely  [1,6]: 


24 


in  which 


A  A 

+  xdx ' ,  or  +  ydy',  for  y=(0;b),  or  x=(0;a), 
respectively,  in  the  case  of  the  rectangular  slit 

d l'  =  i'dl'  = 

A 

+4> ' ad4> 1  in  the  case  of  the  circular  inlet. 


Figure  8.  Scattering  from  the  open  end  and  its  equivalent  problem. 


r  *  position  vector  from  0  to  P  ; 

r'  =  position  vector  from  0  to  mid-point  of  the  rim 

element  dP  {see  Figure  8(b)! 

R  =  [r-r'l  ; 

The  angles  t,  V,  S0  and  0  are  associated  with  the  soft  and  hard  edge 
diffraction  coefficients  Ds  and  0^  which  are  present  in  (31)  and  (32). 
In  particular,  these  edge  diffraction  coefficients  Ds^(ii>,  P;  30, 
in  (31)  and  (32)  are  defined  for  each  point  of  diffraction  on  the  rim 
edge,  and  they  are  associated  with  an  equivalent  edge  of  a  half  plane 
which  is  tangent  to  both  the  rim  edge  as  well  as  the  inlet  surface  at 

each  point  of  diffraction  on  the  rim.  Assuming  that  Q  is  any  such 

point  of  diffraction  on  an  edge  in  an  arbitrary  curved  wedge,  the 
meaning  of  30,  3,  <|>  and  with  respect  to  the  point  of  diffraction  0 
becomes  clear  from  Figures  9(a)  and  9(b)  for  this  wedge  configuration. 
Of  course,  the  walls  of  the  inlets  or  waveguides  are  assumed  to  be 
infinitesimally  thin  in  the  present  work,  so  the  wedge  geometry  in 
Figures  9(a)  and  9(b)  should  be  reduced  to  an  edge  in  a  curved  or  plane 
screen  (n=2  case  in  Figure  9(b))  for  the  circular  and  rectangular 
inlets  in  the  present  instance;  however,  the  meaning  of  <|>,  p,  30  and  3 
is  unchanged  even  when  the  wedge  in  Figures  9(a)  and  9(b)  collapses  to 
an  edge  in  a  screen. 

Finally,  the  equivalent  currents  I6q  and  Meq  are  valid  provided 
one  is  observing  the  fields  of  these  sources  in  directions  away  from 
the  geometrical  optics  incident  and  reflection  shadow  boundaries  [3,6]; 
under  these  conditions  which  are  met  in  practice  (since  directions  of 
incidence  and  scattering  are  restricted  to  60°  or  so  from  the  axis  of 
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the  inlets),  Ds  and  %  for  the  local  half  planes  at  each  point  of 
diffraction  on  the  rim  become  [6]: 

-e"jlT/4 

H  .8)  = - — - 

h  2/ 2-rr k  |/sinBQsinB  | 

(34) 

The  use  of  (31),  (32)  and  (34)  into  (33)  together  with  of  (1;  2(a); 
2(b))  and  =  k1  x  Y^E1  allows  one  to  easily  identify  [S^]  via 
(5)  and  (7).  Thus 


r»M>\ 

sec(— 2“J  +  sec(— yt) 


/3  ■=  Po  0N 

KELLER  CONE 

(a)  (b) 


Figure  9.  Diffraction  by  a  wedge.  The  angles  0?,  8,  V  and  which 
occur  in  the  edge  diffraction  coefficient. 
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Some  typical  results  for  the  fields  scattered  by  only  the  open  end 
of  circular  and  rectangular  inlets  are  shown  in  Figures  10  -  15.  The 
plane  wave  is  incident  on  the  open  end  from  the  direction  (Q1;^),  and 
the  near  field  is  evaluated  either  for  a  fixed  R  and  6r  with  0  <  <j>r  < 
2tt,  or  for  a  fixed  R  and  $r  with  ]9r|  <  60°.  Here,  the  direction  of 
scattering  is  specified  by  the  set  of  angles  ( 0r ; 4> r )  and  R  denotes  the 
distance  from  the  origin  of  the  coordinate  system  (see  Figure  1)  to  the 
near  zone  where  the  scattered  field  is  observed.  Again,  near  zone  or 
near  field  implies  that  the  observation  (or  field)  point  is  located 
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within  the  "radi ati ng"near  zone  of  the  inlet  opening.  It  is  clear  from 

these  figures  that  the  "total"  scattered  field  (ES0  =  xES0  +  yES0  + 

x  y 

zE^°)  generally  contains  a  peak  around  (0r  =  01  and  4>r  =  it  +  4> 1 ) ;  this 
is  to  be  expected  since  it  is  like  a  specular  direction  where  the 
fields  diffracted  from  the  entire  rim  (at  the  open  end)  contribute 
almost  in  phase.  In  these  figures,  all  the  lengths  are  expressed  in 
terms  of  the  wavelength. 

In  order  to  show  the  effects  of  multiple  interactions  of  the  rays 
across  the  aperture,  it  is  useful  to  illustrate  the  on-axis 
backscattered  field  for  a  circular  inlet  as  calculated  previously  in 
[9]  which  includes  up  to  first-order,  second-order,  and  third-order 

diffractions  as  well  as  the  effects  of  all  the  multiple  interactions. 
These  results  as  shown  in  Figure  16  are  calculated  from  an  asymptotic 
high  frequency  ECM  based  approach  and  are  compared  with  an  exact 
solution  based  on  the  Wiener-Hopf  technique.  It  is  seen  that  for 
a/X  >  0.5,  the  backscattered  field  starts  to  settle  down  and  approaches 
(within  2  dB)  the  values  based  on  including  only  first-order 
diffraction.  As  the  electrical  dimensions  of  the  inlet  opening 
increases,  the  results  based  on  using  only  first-order  diffraction  will 
become  more  and  more  accurate.  In  the  computer  code  associated  with 
this  report,  the  first-order  and  the  second  order  edge  diffractions  are 
included  to  show  the  effect  of  the  edges. 


29 


B.  Development  of  [Sgi] 


The  scattering  matrix  [S21]  describes  the  transmission  or  coupling 
of  the  incident  plane  wave  field  into  the  waveguide  modes  as 
illustrated  in  Figure  4  and  described  in  (14)  and  (15). 

For  later  convenience,  let  the  waveguide  in  Figure  4  be  excited  by 
an  electric  current  moment  dpe  which  lies  in  the  exterior  but  near  zone 
region  of  the  inlet  (waveguide)  opening  as  shown  in  Figure  17.  The 
[S21]  will  be  developed  here  for  the  coupling  of  the  fields  of  the 
exterior  source  dpe  into  the  interior  waveguide  region  via  the  inlet 
opening  as  in  Figure  17;  this  more  general  situation  reduces  to  the 
special  case  of  plane  wave  incidence  on  the  inlet  opening  as  in  Figure 
1  if  the  source  dpe  is  allowed  to  receed  to  infinity. 

The  field  E?n  which  is  coupled  into  the  waveguide  region  by  dpe 

in  Figure  17  may  be  expressed  as  in  (8)  by 


E7  =  ll  [A"  e  e 
in  “  1  nm  nm 
mn 


+  i8  z  +iSz 

nm  +  B-  (e  -  ze  )e  nm  ] 
nnr  nm  znm'  J 


(37) 


The  associated  magnetic  field  H7n  is  likewise  given  via  (9)  as 


h:  =  ll  [a;,  (-h'  *  S' 

mn 


in 


nm 


nm 


+j 3  z  _  +je  z 

)  e  nm  +  B"  f-h  )e  nm  ] 
znrrr  nnr  nnr 


(38) 

The  modal  coupling  or  transmission  coefficients  Anm  and  Bnm  in  (37)  and 
(38)  may  be  found  from  the  equivalent  problem  in  Figure  17(c)  which 
illustrates  an  equivalent  surface  and  line  source  distribution  at  z=0 
within  an  Infinite  waveguide  (which  is  an  extension  of  the 
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FIXED  THETA  PATTERN;  NEAR  FIELD  CASE 
EDGE  0  I FFRAC  T I  ON  FOR  CIRCULAR  WAVEGUIDE 
(THI . PHI  I  ==  3Q. 00  0.00 

IRR, THAI  ==  4.00  30.00  A  =  =  1.000 


F I XEO  THETA  PATTERN:  FAR  FIELD  CASE 
EDGE  DETRACTION  FOR  CIRCULAR  WAVEGUIDE 
[THI  .PHI)  ==  30.00  0.00 

R  =  = 


-total  field  due  to  Eg  _ total  field  due  to  E' 


Figure  10.  A  comparison  of  the  far  zone  and  near  zone  edge  diffraction 
from  the  open  end  of  circular  Inlet  for  a  30°  cone  scan. 


DB  PATTERN 


Figure  11.  A  comparison  of  the  far  zone  and  near  zone  edge  diffraction 
from  the  open  end  of  circular  Inlet  for  $  *  0  cut. 


Fixer  THETR  PPITTERN;  NEfiP  FIELD  CH5E 
EOGE'  DIFFRACTION  FOR  CIRCULAR  WAVEGUIDE 
!Thi.r-Hii==  5.00  o.oo 
lRR.THH>=  =  4.00  5.00  R  =  =  1.000 


PHR 

(a) 

FIXEO  THETR  PATTERN;  NEF1R  F1EL0  COSE 
EOGE  0 !  FFRRCT  I  ON  FOR  CIRCULAR  WRVEGU I  DE 
( T  H I , PH  I  1  =  =  5.00  0.00 

(RR. THR)  =s  4.00  30.00  P  =  =  1.000 


PHR 

(b) 

- total  field  due  to  Ej  _ .total  field  due  to  eJ 

Figure  12.  A  comparison  for  0  =  5°  and  8  =  30°  cone  scan  patterns  with 
9'  =  5°  plane  wave  incidence. 
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FIXED  THETR  PRTTERN:  NEAR  FIELD  CASE 
EOGE  DIFFRACTION  FOR  RECTANGULAR  WAVEGUIDE 
[ TH1 . FH 1 1 ==  5.00  0.00 

(RR.  THKt  ==  4.00  30.00  IA,BI=  =  2.000  1.000 


De  PRTTERN 


FIXED  THETR  PRTTERN;  NERR  FIELD  CRSE 
EDGE  DIFFRACTION  FOR  RECTRNGULRR  WAVEGUIDE 
(THI . PHI ) ==  5.00  90.00 

(PR. THR1 ==  4.00  30.00  [R.B)==  2.000  1.000 


Figure  13.  (Continued) 


FIXED  THE  TA  PATTERN:  NEAR  FIELD  CASE 
EDGE  DIFFRACTION  FOR  RECTANGULAR  WAVEGUIDE 
t  T  H I .PHI)  =  =  3U. 00  0. 00 

(RR.THRl==  4.00  30.00  t(A .  B)  =  =  1.770  1.770 


(a) 


FIXED  THETA  PATTERN;  NEAR  FIELD  CASE 
EDGE  0 1 FCRAC  T I  ON  FOR  RECTANGULAR  WAVEGUIDE 
M  H] , PHI ) ==  30.00  0.00 

l HR . T HR)  -  -  4.00  30.00  L A . B)  =  =  2.000  2.000 


(b) 


_ total  field  due  to  Eg  _ total  field  due  to  E^ 

Figure  15.  The  30°  cone  scan  patterns  for  rectangular  inlet  with 

comparable  area  or  comparable  dimension  as  that  of  circular 
inlet  in  Figure  10. 
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FIXED  THEIR  PPTTERN:  FRR  FIELD  COSE 
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Figure  15.  (Continued) 


39 


<0 


*  : 

£ 

So  j 

Ist- ORDER 

UTD 

CO  ™ 
CJ 

oc=fj 

to  : 

0. 

00 

'  '  i  '  ■  '  ■  i  ■  ■  ■  ■  i  ■  ' 

0.25  0.50  0.75 

1.00 

'  1  ■'  ' 
1.25 

Q  /  X 

(a)  first-order 


■n 

l.SO 


a  /X 


(b)  second-order 


Figure  16.  On-axis  backscattered  RCS  (normalized  to  na2)  from  a 
semi-infinite,  hollow,  perfectly-conducting  circular 
cylinder  with  various  terms  included  in  UTD  calculations. 


40 


(c)  third-order 


(d)  all  interaction  terms  included 
Figure  16.  (Continued). 
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(a)  PART  OF[S2l]  DUE  TO  COUPLING  OF  THE 

DIRECT  FIELD  OF  dpTHROUGH  THE  APERTURE 

e 


CORRECTION  FROM  THE  RIM 


(z  =0  PLANE) 

f 


ds 


(c)  EQUIVALENT  PROBLEM 


Figure  17.  Geometry  associated  with  [S21 3  calculation. 
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semi -infinite  inlet  or  waveguide  for  z>0).  The  equivalent  sources  at 
z=0  in  Figure  17(c)  generate  the  same  fields  in  the  waveguide  region 
z<0  as  those  which  exist  in  Figures  17(a)  and  (b)  if  the  equivalent 
sources  are  found  exactly.  Here,  the  equivalent  sources  are  determined 
from  asymptotic  high  frequency  techniques  so  that  the  fields  coupled 
into  the  waveguide  are  approximations  to  the  true  fields  therein. 

These  approximations  are  high  frequency  approximations  which  are 
expected  to  work  well  even  down  to  the  lowest  propagating  mode  in  the 
waveguide.  According  to  the  high  frequency  estimates  based  on  the  GTD, 


(39) 


M 


s 


(40) 


where  and  fr  represent  the  unperturbed  or  geometrical  optics  (GO) 

field  produced  by  dpe  within  the  aperture  region  (in  the  z=0  plane) 

but  in  the  absence  of  the  inlet  structure.  The  additional 

contributions  to  Js  and  must  arise  from  the  diffraction  by  the  edges 

of  the  aperture;  these  edge  diffracted  field  contributions  are  denoted 

by  and  in  (39)  and  (40).  It  is  easily  seen  that  JGG  and  MGG  may 
s  s  s  s 

be  expressed  as: 


-  U  x  Hjp  ] 


e  z=0 


(41) 


(42) 
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where  f e!  ;  h!  )  are  the  (electric;  magnetic)  fields  incident  at  the 

1  dPe  dpe; 

open  end  (z=0)  from  the  external  source  dpe  at  P.  These  incident 


fields  (E^p  ;  )  represent  the  unperturbed  fields  of  dpe  which  exist 

in  the  absence  of  the  inlet.  Clearly  [6,15]: 


'O  i  *  -j  _ 

—  S  X  s  X  dpe 


-jk  ..  e-jks1 
"77  5  x  3pe  — T~ 


The  expressions  in  (43)  and  (44)  are  valid  for  distances  si  which 
correspond  to  dpe  being  in  the  near  zone  of  the  inlet  aperture  (z=0). 
However,  s'*  cannot  be  made  extremely  small  to  where  the  reactive  field 
terms  of  the  type  l/(s^)2  and  1/ ( s^  )3  become  significant;  these  higher 
order  range  dependent  terms  are  ignored  in  (43)  and  (44)  which  pertain 
only  to  the  radiation  fields.  Before  proceeding  to  calculate  the  modal 
transmission  coefficients  Anm,  it  is  convenient  to  decompose  Anpi  as 
follows : 


A"  =  AG0  +  Ad 
nm  nm  nm 


where  AGG  is  the  part  of  A"  due  to  JG0  and  MGG;  whereas,  Ad  is  the 
nm  nm  s  s  nm 

oart  of  A"  which  is  produced  by  Jd  and  Md.  It  is  now  an  easy  matter 
nm  s  s 

to  find  AGG  from  JGG  and  by  employing  the  results  of  Appendix  I 
which  indicate  the  manner  in  which  electric  and  magnetic  current 
sources  excite  modes  inside  a  waveguide  [13].  Thus,  from  (A-8)  of 
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Appendix  I,  the  GO  part  of  the  TE  transmission  (or  coupling) 

coefficients  A  are  given  by: 

nm  3  J 


-II  ds[(e  e 
(aperture(z=o))  nm 


•J'B  z 


-  I  ^  I 


)*(r)-(h  +zh  )e 
'  '  s  ’  K  nm  znm; 


J^nmz. ruGO 


(MsbU)] 


2  I  I  ds  z  •  (e  x  h  ) 
(aperture(z=0))  nm  nm 


In  a  similar  fashion,  one  may  express  the  TMnm  transmission  (or 
coupling)  coefficients  B”  by 

■  *■'  '  n  m 


b"  =  BG0  +  Bd 
nm  nm  nm 


(47) 


where  B  u  is  obtained  via  (A-8)  as 
nm 


rn  -II  dsK«  +«  )e’J@nn,Z-C  JG0)-Ch  e”J  B"mZ)  -  C  MG0)  ] 

B“  =  (aperture(z-o))  nm  zm  5  nm  5 


2  J  /  ds  z  •  (e  x  h  ) 
(aperture(z=0) )  nm  nm 


It  is  not  easy,  in  general  to  calculate  Ad  and  Bd  from  Jd  and  Md  in 
J  3  nm  nm  s  s 

the  aperture  (z=0);  however,  Ad  and  Bd  may  be  calculated  indirectly 
r  nm  nm 

via  an  approach  based  on  the  ECM  which  employs  equivalent  magnetic  line 
and  line  dipole  sources  on  the  rim  of  the  aperture  (z=0).  These 

o  <*  j)  H'' 

equivalent  sources  of  strengths  M  i '  and  M  z  corresponding  to  the 
magnetic  line  and  line  dipole  rim  currents  radiate  "inside  the  infinite 
waveguide"  (corresponding  to  the  inlet  and  it's  extension  for  z>0)  as 
in  Figure  17(c).  Since  the  equivalent  rim  currents  must  lie  on  the 
"waveguide  walls"  (see  Figure  17)  in  this  case,  only  the  equivalent 
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magnetic  type  rim  currents  are  thus  required  as  any  equivalent  electric 
type  rim  currents  which  are  tangential  to  the  waveguide  walls  would  be 
shorted  out  because  the  walls  are  perfectly-conducting.  This  situation 


is  in  contrast  to  the  exterior  radiation  or  [S^]  calculation  via  the 

equivalent  electric  and  magnetic  rim  currents  I  and  M  which  radiate 

eq  eq 

in  "free  space"  as  shown  in  Figure  8(b).  Also,  I  and  M  of  (31)  and 

eq  eq 

(32)  contain  Keller's  6TD  edge  diffraction  coefficients  D$  and  ; 

whereas,  for  reasons  which  are  indicated  below,  the  equivalent  rim 
l  £d 

currents  M  and  M  are  defined  in  terms  of  the  Ufimstev  type  edge 
diffraction  coefficients  and  0^,  respectively,  as  follows*  [61: 
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(50) 


where  DU  and  Dd  are  given  by  [6]: 
s  n 


The  fields  and  '  are  multiplied  by  a  factor  of  1/2  in  (49) 

and  (50)  to  account  for  the  presence  of  the  waveguide  walls  on  which  nf 

j) 

and  M  are  located:  this  aspect  is  discussed  in  [6], 
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(51) 

It  is  necessary  to  use  the  Ufimstev  diffraction  coefficients  0^  and  ojj, 

respectively,  in  (49)  and  (50)  instead  of  the  Keller  edge  diffraction 

coefficients  0  and  0.  due  to  the  fact  that  the  truncation  of  the 

s  h  s 

“Pn 

and  M“  sources  by  the  rim  or  boundary  (end  points)  of  the  surface 

integrals  in  (46)  and  (48)  over  the  inlet  aperture  (z=0)  gives  rise  to 

some  contributions  to  the  edge  diffracted  fields  which  propagate  into  the 

waveguide.  Hence,  the  remaining  contributions  to  the  edge  diffracted 

field  must  come  from  h  according  to  the  PTD  ansatz  instead  of  from 

D  ,  which  yields  the  total  first  order  diffracted  field  [6],  Clearly, 
s  j  n 

the  use  of  D$  ^  instead  of  ^  in  (51)  would  erroneously  include  some 
diffraction  effects  twice  when  the  effects  of  (49)  and  (50)  are 
combined  with  (46)  and  (48).  In  (49),  (50)  and  (51),  the  incident  ray 

a 

from  dpe  makes  an  azimuthal  angle  about  the  unit  edge  tangent  at 

any  point  of  diffraction  on  the  rim  (see  Figures  9(a)  and  9(b)). 

Likewise,  in  ( 49 ; 50 ; 51 )  corresponds  to  an  azimuthal  angle  which  the 
J 

ray  diffracted  from  the  edge  along  the  s^  direction  into  the  waveguide 

A 

region  makes  about  l'  at  the  same  point  of  diffraction  (note:  3nmj.  and 

<j/)m  are  like  3  and  in  Figures  9(a)  and  9(b);  likewise  3  and  i|»' 

J  O 
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also  have  the  same  meaning  as  in  those  figures.)  It  is  important  to 
also  note  that  the  direction  of  the  ray  diffracted  into  the  waveguide 
must  be  chosen  to  correspond  to  the  "modal  ray"  direction  [6].  The 
circular  waveguide  modes  can  be  asymptotically  approximated  in  the 
circular  inlet  case  by  a  set  of  axially  divergent  and  convergent 
conical  modal  ray  systems  sufficiently  far  from  the  waveguide  axis.  If 
one  considers  the  reciprocal  problem  of  the  radiation  into  the  exterior 
region  by  the  modes  incident  at  the  open  end  from  within  the  waveguide 
region  (z<0),  then  it  is  only  the  upgoing  or  axially  divergent  (and  not 
axially  convergent)  conical  modal  ray  system  which  constitutes  the 
field  "incident"  on  the  edge;  this  incident  field  diffracts  from  the 
rim  edge  (at  z=0)  to  contribute  to  the  exterior  field.  The  of 

J 

(49; 50)  in  the  original  problem  is  associated  with  the  TE  and  TM 
conical  modal  rays  which  are  launched  via  edge  diffraction  and  which 
exactly  correspond  to  the  axially  divergent  conical  modal  rays  in  the 
reciprocal  problem  except  that  these  rays  are  now  reversed  with  respect 
to  those  in  the  reciprocal  problem.  Alternately  phrased,  \|/!m  of  (49;50) 

\J 

in  the  original  problem  is  associated  with  the  axially  convergent  conical 

model  rays  which  are  launched  into  the  waveguide  region  due  to  edge 

diffraction  effects  contained  in  and  M^.  Since  one  deals  with  only 

one  set  of  conical  modal  rays  (i.e.,  of  the  axially  convergent  type  in 

this  situation)  the  sums  in  (49)  and  (50)  reduce  to  a  single  term  for 

each  TE  and  TM  case,  respectively;  i.e.,  N=1  in  (49)  and  (50)  for  the 
nm  nm 

circular  inlet  configuration.  On  the  other  hand,  the  rectangular 
waveguide  modes  can  be  expressed  as  a  set  of  four  modal  plane  wave  or  ray 


48 


fields  everywhere  inside  the  waveguide  including  the  edges  at  the  open 

end.  These  rays  also  contain  a  pair  of  rays  that  propagate  toward  each 

of  the  four  edges  as  well  as  a  pair  of  rays  that  propagate  away  from 

those  edges.  In  the  present  1  calculation  ,  it  is  of  interest  to 

associate  in  the  case  of  the  rectangula"  inlet,  to  rays  launched 

into  the  waveguide  region  ( z <0 )  via  edge  diffraction  of  the  field 

incident  from  dpe  (which  lies  in  the  exterior  region  as  in  Figure  9). 

Thus,  is  the  azimuthal  angle  (about  l')  associated  with  the  edge 

diffracted  rays  which  are  identified  with  a  pair-  of  rectangular  waveguide 

modal  rays  that  propagate  away  from  the  edges.  Clearly,  since  there 

exists  a  pair  of  such  modal  rays  in  this  [S21l  calculation  for  the 

rectangular  inlet,  the  summing  index  in  (49)  and  (50)  goes  from  j=l  to 

j=2;  i.e.,  N=2  in  (49)  and  (50)  for  the  rectangular  inlet  case.  The 
a  id 

radiation  of  M  and  M  0f  (49)  aruj  (50)  into  the  waveguide  region  is 
given  as  before  via  (A-8)  of  Appendix  I.  Thus,  the  coefficients  A^m  and 

A  g  O/j 

r  of  the  TE  and  TM  modes  excited  by  M  and  M  are  explicitly  given 
unm  nm  nm  J  3 

by 


l  f  (d£ '  [ h  •  +  zh'  *M*dz]e  J'  nm  ) 

j=l(rim(z=0))  nm  J  znm  j  J 


nm 


2  //  ds  z 
(aperture  (z=0)) 


(i  *  R  ) 

'  nm  nm 


(52) 


and. 


N 

l  f  (dt 1  [h  « 
j=l(rim(z=0))  ™ 


nti']e‘j8nn’Z 

J  J 


) 


2  //  ds  z 
(aperture  (z=0)) 


fi  ) 

nm7 


(53) 
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C.  Near  Zone  Coupling  Matrix  [S21] 


The  scattering  matrix  [S21 ]  can  easily  be  identified  via  (14),  (15), 
(45),  (46),  (47),  (48),  (52)  and  (53).  It  is  noted  that  M®°,  3®°,  M* 

^  ^  J 

and  contain  information  on  A„  and  A,  via  E1  and  that  are  present 
j  0  <P 

in  these  quantities  (also  see  (1)  and  (2)).  Recall,  of  course,  that 
the  present  development  can  also  deal  with  a  near  source  dpe  type 
illumination  as  indicated  in  (43)  and  (44)  in  which  case  the  coupling 
matrix  [S21]  must  be  rePlace<1  by  ^S2i^  defined  via  the 

relationship 


[A'  ] 

L  nmJ 

[s'  ]  [s'  ]  [s'  ] 

L  xnmJ  L  ynmJ  L  znmJ 

dp 

dpex 

[B“  ] 

L  nmJ 

[s  ]  [s  ]  [S  ] 

L  xnmJ  L  ynmJ  L  znmJ 

dpe* 

Kez 

where 


[S21'1  = 


[s  ]  [s  ]  [s  ] 

L  xnmJ  L  ynmJ  L  znmJ 

[s  ]  [s  ]  [s  ] 

L  xnmJ  L  ynmJ  L  znmJ 


and 


(54) 


(55) 


—  — 

x  •  dp 

dp 

c 

A 

dp®* 

= 

y  *  dp 

dP^ 

e 

A 

ez 

z  •  dp 

e 

If  the  source  dpe  receedes  to  infinity,  then  corresponding  to  a 
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locally  plane  wave  type  of  illumination.  Also,  [S^]  [dpg]  on  the  right 
side  of  (54)  should  approach  [S21]  EE1]  which  occurred  on  the  right  side 
of  (14)  when  dpg  receedes  to  infinity.  Clearly  then. 


[A'  ] 

L  nmJ 

=  [sL] 

dp 

dpf 

-4  [S_-  ] 

Ae 

A® 

[B‘  1 

L  nmJ 

21 

dpey 

*ez 

21 

4> 

if  dp  receedes 
e 

to  infinity 

(56) 


Some  numerical  results  for  [S21]  will  be  indicated  in  terms  of 

<v 

[S  2]  since  [S12]  is  obtained  from  [S21]  via  reciprocity. 


D.  Development  of  [S^]  Using  Reciprocity 

The  scattering  matrix  [S12J  describes  the  manner  in  which  the 

modal  fields  propagating  within  the  waveguide  region  (z<0)  radiate  into 

the  exterior  region  from  the  open  front  end  (at  z=0)  as  in  Figure  5. 

The  observation  point  can  be  in  the  near  zone  of  the  inlet  aperture  as 

mentioned  previously.  This  matrix  [S12]  was  defined  earlier  in  (10) 

ro 

and  (12).  Actually  the  problem  of  determining  EXfy>z  of  (10)  in  the 
near  zone  of  the  inlet  aperture  can  be  seen  to  be  related  to  the 
problem  of  determining  Anm  and  Bnm  due  to  the  source  dpe  at  the  point  P 
in  the  near  zone  as  In  Figure  17  via  the  reciprocity  theorem  for 
electromagnetic  fields.  Thus,  [S12]  of  (12)  may  be  found  from  [S21]  of 
(55)  via  the  reciprocity  principle  as  illustrated  below. 

Let  an  electric  current  moment  dpe  at  P  in  the  near  zone  of  the 
inlet  aperture  generate  the  fields  (E  ;  H  )  everywhere1-  in  the  presence 
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of  the  inlet  structure  which  is  asumed  to  be  a  semi -infinite  hollow 


waveguide  as  shown  in  Figure  18.  Also  let  a  source-free  waveguide 
modal  field  E*  propagate  from  the  region  z<0  to  the  inlet  opening  at 

z=0,  and  let  this  mode  generate  the  total  fields  (E  ;  H  )  after  it 

-po  “po 

scatters  from  the  open  end.  Furthermore,  let  (E  ;  H  )  exist 
everywhere*  with  dpe  absent.  Now,  an  application  of  the  reciprocity 

"A  “p  -p0  “PO 

theorem  to  the  pair  of  fields  (E  ;  H  )  and  (E  ;  H  )  within  the  region 
V  bounded  by  the  surfaces  S,  Sq  and  the  surface  Z  at  infinity  (see 
Figure  18)  yields: 

If  [Erox  He  -  Ee  x  Hro]  •  n  ds  =  dpe  •  Er0(P)  .  (57) 

S+S0+E 

Next,  employing  the  boundary  conditions  [nxE  ]$  =  0,  [nxE  =  0  on 
the  surface  which  encapsulates  the  waveguide  walls,  and  also  enforcing 
the  radiation  conditions  on  (Er0;Hr0)  and  (EP;He)  over  the  surface  Z  at 
Infinity  in  (57)  leads  to 

-  II  He-Mpqds  +  //  Ee.j®qds  =  dpe  •  Ero(P)  ,  (58) 

S0  S0 

where  Meq  and  Jeq  constitute  the  field  quantities  associated  with  Ero 
s  s 

and  Hro  Inside  the  waveguide.  In  particular 


*Here,  "everywhere"  includes  the  regions  which  are  both  interior  and 
exterior  to  the  inlet. 
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CO 


Figure  18(a).  Geometry  for  an  application  of  the  reciprocity  theorem 
which  relates  [S^l  and  CSgi  1  • 
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H  [Erox  n]: 


i  [n  x  Hro], 


Furthermore, 


=  (E+  +  l  R  E")  A+  , 
1 P  q  PUJ  P 


-  (H+  +  I  R  A+ 

v  p  q  pq  p 


The  fields  (E+;  H+)  are  as  defined  in  (A-2)  of  Appendix  I.  Here 
P  P 

(E*;  Hp)  Is  the  pth  mode  Incident  at  the  open  end  of  the  waveguide  as 

in  Figure  5.  The  single  subscript  p  refers  to  a  double  mode  index  n,m 

as  mentioned  in  Appendix  I;  likewise,  the  subscript  q  refers  to  another 

pair  of  modal  Indices  1,j,  for  example,  where  1  and  j  can  be  different 

from  n  and  m.  The  Incident  modal  fields  (E*;  H*)  generate  the 

P  P 

reflected  fields  (E  R  E~  ;  £  R  H")  after  scattering  at  the  open 

q  qp  q  q  qp  q 

end,  and  Rqp  denotes  the  reflection  coefficient  or  the  amplitude  of  the 
modal  field  which  is  reflected  back  into  the  waveguide  after  a  pth 
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mode  is  incident  at  the  open  end.  One  can  also  write  Ee 
as  follows: 


and  He 

So 


z  A„Er 

i  1  % 


(63) 


=  Z 


a'h“ 

l  l 


(64) 


where  the  subscript  i  once  again  refers  to  a  double  mode  index  rs,  for 
example,  as  in  Appendix  I.  Incorporating  (59)-(64)  into  (58),  and 
making  use  of  the  orthogonality  property  of  the  waveguide  modes  gives: 


dp  •  Ero(P)  =  -  2Af  A-  //  E  x  H 
Ke  p  p  y  pt  pt 

^o 


zds 


(65) 


In  matrix  notation,  it  is  clear  that  (65)  can  be  expressed  as: 


Er°(p)  Ero(p}  ErO(P) 
x  y  z 


dpex 

dPey 

dpez 
_  _i 


([A;])T(-2[//E"pt  »  Hpt-Sds])([A-]). 
0 

(66) 


It  is  noted  that  Er0  contains  A^  since  Er0(P)  is  produced  by  the 

radiation  of  A+E+.  Furthermore,  it  is  obvious  that  (10)  and  (54)  can 
P  P 

be  expressed  more  compactly  as: 
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e™(p) 

Ero(P) 

y 

Ero(P) 

z 


ts12] 


(67) 


and 


—  — 

dp 

[A“] 

L  PJ 

=  [S21] 

dpex 

dpey 

Kez 

where 


(68) 


Incorporating  (67)  and  (68)  into  (66)  yields: 


(69) 


T 

cs12f 


dp 

dpex 

dpe* 

ez 


—  — 

T 

dp 

[A+] 

L  PJ 

(-2[//  Ept  «  Spt  •  i«])[Sa] 

dpex 

dp** 

—  — 

0 

ez 

(70) 


The  relationship  between  [S 1 2 H  and  C$21]  now  becomes  evident  from  (70); 
namely: 
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(71) 


[S12f  *  (-*[//  Epl  x  H‘pt  •  ids])  [S21l  . 

O 

rv  <V 

Therefore,  one  can  find  [S12]  from  a  knowledge  of  [S21V,  fIS2l  1  was 
obtained  previously  in  part  B  of  this  section. 

Typical  near  zone  radiation  patterns  of  open-ended  rectangular  and 
circular  waveguides  which  are  calculated  by  the  above  procedure  are 
indicated  in  Figures  19  to  33.  If  the  near  zone  field  point  is  allowed 
to  recede  to  the  far  zone  then  the  above  procedure  also  furnishes  the 
far  zone  radiation  patterns.  Indeed,  the  latter  are  shown  in  Figures  19 
to  23  for  the  open  ended  circular  waveguide  problem;  these  principal 
plane  far  zone  radiation  pattern  calculations  are  compared  with  those 
obtained  via  the  rigorous  Wiener-Hopf  method  [161  in  Figures  19  to  23; 
the  comparison  is  seen  to  be  quite  good  in  these  figures.  It  is  noted 
that  while  the  Wiener-Hopf  method  provides  a  useful  check  on  our 
solutions,  the  present  procedure  is  far  simpler  than  that  based  on  the 
Wiener-Hopf  method.  There  is  presently  no  corresponding  rigorous 
analytical  solution  available  in  the  literature  for  the  far  zone 
radiation  pattern  of  an  open-ended  rectangular  waveguide  as  there  is  for 
the  circular  waveguide  case;  hence,  the  rectangular  waveguide  far  zone 
radiation  patterns  are  shown  without  any  comparisons.  Explicit 
expressions  for  the  far  zone  radiation  patterns  for  the  circular  and 
rectangular  inlets  are  presented  for  convenience  in  Appendix  III. 
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.present  solution  x  x  x  x  Weiner-Hopf  [16] 


Figure  19.  E*  radiation  pattern  due  to  a  TEqi  mode  In  an  open-ended 
circular  waveguide. 
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corresponding  to  ^9  far  zone  patterns 


Figure  24.  The  near  zone  radiation  patterns  of  TE^  modes 
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Figure  29.  The  near  zone  radiation  pattern  of  TM02  mode 
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Figure  30.  Near  zone  radiation  pattern  of  an  open-ended  rectangular 
waveguide.  kaa12.58,  aa2b,  TEjq  mode,  ^p»0°,  kR»62.90 


Figure  31.  Near  zone  radiation  pattern  of  an  open-ended  rectangular 
waveguide.  ka=12.58,  a=2b,  TE^q  mode,  <j>r=90°,  kR=62„90. 


Figure  32.  Far  zone  radiation  pattern  of  an  open-ended  rectangular 
waveguide.  ka=12.58,  a=2b,  TE^q  mode,  <t>r=0° . 
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Figure  33.  Far  zi 
wavegi 


E.  Development  of  CS22I 

As  indicated  earlier  in  Section  I,  represents  a  modal 

reflection  coefficient  matrix  which  is  associated  with  the  interaction 

illustrated  in  Figures  5  and  34.  In  particular,  is  made  up  of 

elements  which  correspond  to  the  modal  amplitudes  (which  are 

proportional  to  the  reflection  coefficients)  of  the  TE  and  TM  modes 

nm  nm 

that  are  reflected  back  into  the  waveguide  region  (z<0)  when  either  a 

TE  or  a  TM  mode  is  incident  on  that  open  end  from  the  region  z<0. 

pq  pq 

Once  again,  GTO  and  ECM  will  be  employed  to  calculate  these  modal 
reflection  coefficients.  It  was  mentioned  in  part  B  (of  this  section) 
dealing  with  the  development  of  1  that  the  modal  fields  inside  a 
circular  or  rectangular  waveguide  could  be  decomposed  into  fields  which 
are  associated  with  a  set  of  "modal  rays".  There  exists  only  one  set  of 
rays,  namely  the  axially  divergent  modal  rays,  which  are  "incident"  on 
the  circular  inlet  edge;  these  incident  rays  can  partly  diffract  from 
the  edge  to  propagate  back  into  the  waveguide  region  as  a  set  of  axially 
convergent  modal  rays  to  form  the  reflected  modes.  On  the  other  hand, 
there  is  a  pair  of  such  modal  rays  which  are  incident  on  the  edge  of  a 
rectangular  inlet;  again,  these  rays  can  diffract  at  the  edge  to  enter 
back  into  the  waveguide  as  a  pair  of  rays  which  exit  or  leave  the  edge 
to  form  the  reflected  modes.  It  is  noted  that  the  pair  of  rays  in  the 
rectangular  inlet  case  contain  oppositely  traveling  plane  wave 
components  with  respect  to  each  edge;  consequently  these  oppositely 
traveling  components  flip  directions  upon  undergoing  diffraction  at  the 
edges  to  convert  the  incident  pair  of  modal  rays  into  the  reflected  pair 
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(a)  GEOMETRY  FOR  THE  REFLECTION  COEFFICIENT 
[S22]  MATRIX 


( Z  *  0  PLANE) 

(b)  EQUIVALENT  PROBLEM 


Figure  34.  Geometry  for  the  reflection  coefficient  matrix  [S22] 
calculation. 
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of  modal  rays.  One  may  define  a  set  of  equivalent  magnetic  rim  currents 

Q  ~  0  {]''  f  )n 

of  strengths  1  and  MR  z  which  produce  the  r  modal  reflected  fields 
(E~;  H~)  when  a  set  of  s'-*"1  modal  fields  (E*;  H*)  is  incident  at  the  open 
end  of  the  inlet.  The  fields  entering  back  into  the  waveguide  as 


reflected  modes  are  produced  entirely  via  edge  diffraciton  of  (E+;  H+) 
as  the  ray  picture  reveals;  hence,  the  equivalent  magnetic  line  and 
line  dipole  rim  currents  and  can  be  expressed  as  in  (49)  and  (bO) 
but  with  the  Ufimstev  edge  diffraction  coefficients  in  those  equations 
replaced  by  the  complete  (Keller  type)  edge  diffraction  coefficients 
and  D  ,  respectively,  as  follows: 


N 

r 

/,  "R. 
J=1  J 


m: 


n 

l 

j=i 


Wi-'V-W  W  («s*J 
-zn  V  Jir !  "7  : - ? 


!/ 


si n  0S  sin3p 

j  j 


(72) 


and 


1  Ds(»j-l|,j;  V  Brj) 

jk  |s^*z|  |/^  1  1 


=  l  ■  )'  ( 

R  J-l  J  J-l 


sm3$  sin3p 

j  J 


(73) 


where  0  , 

s,h 


in  the  above  equations  are  given  (as  in  (34))  by: 


-jtr/4 


h 


?/iiT k  |  /si  nB  s  i  n  3  p 


sec 


r  s 

2 


r  s 

lb  .  +  0  ■ 

J 


+  sec  [ 


(74) 
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Following  the  development  of  part  B,  the  currents  and  M^d  generate 
the  reflected  modal  amplitude  associated  with  the  reflected  modal 
field  (i”;  H“)  as  in  (52)  and  (53);  thus: 


l  I  (d£'[h  •  Mp.Jt*  +  zh  •  M^zl  e 

j=l  (rln(z-O))  m  RJ 


2  //  ds  z  •  (e  x  h  ) 
/  *  ,  rtvv  nm  nm 
(aperture  (z=0)) 


-j&z 


I  /  (dt'[h  •  M&.i'  ]  e‘J6™Z) 

J-l  (rlmh-0))  1  nm  «3  J  1 

2  JJ  ds  z  •  (e  *  h  ) 
nm  nm 


where 


[  A+  ] 

1  pqJ 


[B+  ] 

L  pqj 


(77) 


(78) 


0  OA  4. 

It  is  clear  that  and  Mp^  in  (75)  and  (76)  contain  A$  which  are  the 
modal  amplitudes  of  (i*;  H*)  via  (72)  and  (73);  hence,  the  elements  of 
the  reflection  coefficient  matrix  [S^]  may  be  easily  obtained  (even  by 
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inspection)  from  (16),  (17),  (75)  and  (76)  together  with  (73),  (74), 

(77)  and  (78). 

The  modal  reflection  coefficients  or  the  elements  of  [S22]  for  the 
circular  inlet  configuration,  which  are  calculated  by  the  equivalent 
current  method  (ECM)  described  above  are  compared  with  those  based  on 
the  exact  Wiener-Hopf  solution  in  Figures  35  -  38.  This  comparison  in 
Figures  35  -  38  is  also  quite  good.  Once  again,  it  is  noted  that  the 
exact  Wiener-Hopf  solution  [17]  is  far  more  complicated  than  the  present 
one  based  on  the  ECM.  It  is  noted  that  the  effects  of  the  multiple 
diffractions  of  the  rays  across  the  aperture  are  included  in  the 
evaluation  of  the  equivalent  currents  for  the  ECM  analysis  of  [S22] 
pertaining  to  the  circular  inlet  cases  as  shown  in  Figures  35  -  38.  The 
details  of  the  procedure  for  including  these  multiple  interactions  can 
be  found  in  [9];  that  approach  has  been  used  here  as  well  with  the  final 
expression  given  in  Appendix  IV.  In  Figures  39  -  40,  results  for  the 
modal  reflection  coefficients  are  shown  for  open  ended  rectangular 
waveguides  (or  inlets),  and  a  comparison  of  one  such  calculation  with 
measurement  [181  is  shown  in  Figure  41  for  the  dominant  TEiq  rectangular 
waveguide  mode.  The  comparison  with  measurements  in  Figure  41  is 
reasonable;  one  would  expect  this  comparison  to  improve  as  the  frequency 
increases.  The  measurements  in  Figure  41  were  performed  at  a  frequency 
for  which  the  waveguide  dimensions  are  quite  small  in  terms  of  the 
wavelength,  thereby  straining  the  high  frequency  ECM  solution;  it  is 
possible  that  multiple  interactions  of  waves  across  the  aperture  can 
become  important  at  such  lower  frequencies;  on  the  other  hand,  the  ECM 
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Figure  41.  Comparisons  of  the  calculated  and  measured  reflection 

coefficients  in  an  x-band  and  ku-band  open-ended  rectangular 
waveguide. 
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solution  can  become  more  cumbersome  if  these  multiple  wave  interactions 
are  included  for  the  rectangular  inlet  case.  The  effect  of  multiple 
wave  interactions  across  the  aperture  is  more  easily  calculated  in  the 
case  of  the  circular  inlet  geometry. 


F.  Development  of  [Sr 1 

The  matrix  [Sr ]  contains  elements  which  are  the  reflection 

coefficients  of  the  modes  reflected  from  the  termination  at  z=-L  in 

Figure  6  when  a  set  of  modes  is  incident  upon  this  termination.  Let  an 

rth  modal  field  (E~;  H“)  of  the  TE^  or  TM^  type  he  incident  on  the 

termination  at  z=-L  in  Figure  6,  and  let  (E*;  H*)  denote  the  fields  of 

an  sth  TE  or  TM  mode  reflected  from  that  termination.  The  matrix 
s  s 

[S*]  then  relates  the  fields  (E  ;  H*)  reflected  from  ther  termination  at 
r  s  s 

z=-L  to  the  incident  fields  (L“;  H"). 

For  an  impedance  termination,  the  matrix  elements  of  [Sr ]  can  be 
trivially  obtained  by  enforcing  the  impedance  boundary  condition  (at 
z=-L)  on  the  total  rields;  namely: 


[E  -  zz  •  E  =  Z  z  x  H]  (79) 

s  z=-L 

with 


«;  *  s  »;  s 

s 

(80) 

■+  + 

Ar  +  l  \  As 
s 

(81) 
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modes 


in  which  z  (E+;  H+)  are  the  reflected  fields  of  all  the  "s" 
s  Sf h  S 

when  only  an  r  set  of  modal  fields  (E^;  H’)  is  incident  on  the 
impedance  termination  at  z=-L.  Here,  Z$  denotes  the  value  of  the 

surface  impedance.  The  elements  of  [Sp]  as  defined  in  (18)  and  (19) 
can  be  identified  by  simply  incorporating  (80)  and  (81)  into  (79),  and 
by  invoking  the  orthogonality  property  of  the  waveguide  modes.  It  can 
be  easily  verified  that  [Sp]  is  a  diagonal  matrix. 

For  a  planar  isotropic  homogeneous  dielectric/ferrite  at  Z=-L 
which  is  characterized  by  the  permeability  urp0  and  permittivity  ere0, 
the  elements  of  the  matrix  [Sr]  can  be  found  as  follows: 

For  TM  modes: 

0  -  e  0 

„  nm  r  nm 
see  = - 

1  0  +  e  0 

nm  r  nm 


For  TE  modes: 


w  3  -  0 

r  nm  nm 


u  0  +  3 

r  nm  nm 


where 


and 


s™  '  V  koVr  -  kc 
■  V  ko  '  kc 
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for  a  rectangular  waveguide. 


for  circular  waveguide 
TM  modes. 


for  circular  waveguide 
TE  modes. 


Eor  many  materials  of  interest,  the  relative  permeability  ur  «  1, 
therefore  the  reflection  from  the  termination  for  TE  mode  can  be 
simplified  as 


S 


hh 

r 


3  -  3 

nm  nm 


3  +3 

nm  nm 


The  reflection  coefficients  for  blade/disk  and  hub  structure 
terminations  in  circular  inlets  are  slightly  more  complicated  and  are 
presented  in  Appendix  VI. 


IV.  NUMERICAL  RESULTS  AND  COMPARISONS  WITH  MEASUREMENTS  FOR  THE 
SCATTERING  BY  THE  COMPLETE  INLET  CONFIGURATIONS  OF  FIGURE  1. 

In  the  previous  section,  the  various  generalized  scattering 
matrices  were  developed  corresponding  to  the  canonical  scattering  events 
depicted  in  Figures  3-6;  also,  numerical  results  illustrating  the 
behavior  of  these  scattering  mechanisms  were  included.  Comparisons  with 
exact  calculations  or  measurements  were  made  whenever  reasonably 
possible  to  indicate  the  accuracy  of  the  mathematical  expressions  which 
were  developed  to  describe  these  scattering  mechanisms;  those 
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expressions  were  derived  via  asymptotic  high  frequency  techniques  (such 
as  GTD,  UTD,  ECM  or  modified  PTD).  Asymptotic  high  frequency  techniques 
by  their  very  nature  are  approximate;  they  are  expected  to  become 
increasingly  more  accurate  with  increase  in  frequency.  However,  most 
asymptotic  high  frequency  solutions  are  known  to  work  reasonably 
accurately  almost  down  to  frequencies  where  the  characteristic  dimen¬ 
sions  of  the  radiating  object  are  not  much  larger  than  a  wavelength; 
this  is  also  observed  in  the  present  study.  In  this  development  it  is 
necessary  that  the  inlet  cross -section  be  large  enough  to  where  at  least 
the  dominant  mode  can  be  excited  within  the  inlet  waveguide  region. 

Also,  because  of  the  complexity  of  obtaining  the  eigenvalues  in  a 
circular  inlet,  the  computer  program  for  obtaining  numerical  results  for 
the  scattered  fields  can  now  handle  a  circular  inlet  with  radius  up  to  5 
wavelengths,  or  include  up  to  240  modes  (with  the  understanding  that 
there  are  two  independent  azimuthal  wave  functions  cosn<|>  and  s i nn<j>  in 
each  waveguide  mode).  For  rectangular  inlets,  there  is  no  such 
restriction,  since  the  calculation  of  the  eigenvalues  is  very  simple. 
However,  the  calculation  of  scattering  from  extremely  large  inlets  is 
still  not  practical  due  to  the  fact  that  numerous  modes  have  to  be 
included  in  the  calculation.  Further  studies  in  selecting  significant 
modes  and  different  useful  high  frequency  representations  of  waveguide 
modes  are  necessary  in  order  to  get  a  more  efficient  way  to  estimate  the 
scattering  of  large  waveguides  [19,20]. 

The  calculated  numerical  results  shown  In  Figures  42  -  48  include 
the  effect  of  rim  scattering  (i.e.,  the  solid  line  - ),  the  field 
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scattered  by  the  termination  (i.e.,  the  long  dashed  line _ )  and 

the  total  scattered  field  from  the  inlet  which  is  a  vector  sum  of  the 
rim  scattered  field  and  the  scattered  field  from  the  termination  (i.e., 

the  short  dashed  line _ ).  In  these  examples,  the  terminations  are 

a  planar  conducting  surface  and  a  dielectric  termination.  The  short 
circuited  termination  is  a  special  case  of  the  impedance  termination 
with  the  value  of  the  surface  impedance  being  zero. 

For  the  effect  of  rim  scattering,  the  scattered  field  peaks  up  in 
the  neighborhood  of  the  specular  region  (i.e.,  close  to  9r  =  e1, 

<j>r  =  if*1  +ir )  as  is  to  be  expected  from  the  GTD  point  of  view.  Rut  for  the 
field  scattered  by  the  termination,  there  are  no  special  rules  as  to 
where  the  total  field  should  peak  up.  The  field  scattered  by  the 
termination  comes  from  a  combination  of  the  radiation  from  different 
modes,  and  the  phase  between  the  modes  will  change  the  radiation  pattern 
dramatically.  However,  if  one  is  only  interested  in  a  specific  region, 
one  can  pick  up  the  modes  which  contribute  most  significantly  to  that 
region  using  the  concepts  in  [19,20].  In  the  latter  instance,  the 
distribution  of  the  field  within  the  inlet,  which  radiates  out  to 
correspond  to  the  field  scattered  by  the  termination  is  controlled  by 
the  interaction  of  the  propagating  modes.  Therefore,  when  the  length  of 
the  inlet  termination  is  changed,  the  energy  will  be  redistributed 
within  the  various  modes.  For  the  case  of  a  short  circuited 
termination,  the  scattered  field  from  the  termination  is  usually  greater 
than  the  field  scattered  from  the  rim  except  in  the  specular  region.  Ry 
changing  the  impedance  of  the  termination,  from  zero  (short  circuit 
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Figure  42.  Contributions  to  the  4>a180°  patterns  for  circular  inlet. 
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Figure  43.  Contributions  to  the  4>=0°  patterns  for  circular  inlet. 
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Figure  44.  Contributions  to  the  30°  cone  scan  patterns  for  circular 
Inlet. 
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Figure  45.  Contributions  to  the  30°  cone  scan  patterns  for  circular 
inlet  with  almost  matched  termination.  (er  =  1-0. lj) 
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Figure  46.  Contributions  to  the  30°  cone  scan  patterns  for  rectangular 
inlet  with  size  (A ,B )  =  (2,1). 
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Figure  47.  Contributions  to  the  30°  cone  scan  patterns  for  a 
rectangular  inlet  with  size  (A,B)  =  (1,2). 
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(a)  On-axis  RCS  (normalized  ira2 )  of  a  hollow,  perfectly 
conducting  finite  length  circlar  open  at  the  front 
end  and  closed  (shorted)  at  the  hack  end. 

Figure  49.  On-axis  RCS. 
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(b)  On-axis  RCS  (normalized  wa2)  of  a  hollow,  perfectly 
conducting  finite  length  circular  cylinder  open  at 
both  ends. 

Figure  49.  Continued. 
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case)  to  one  which  has  a  non-zero  real  part  (corresponding  to  a  lossy 
impedance),  one  can  drop  the  level  of  the  field  scattered  by  the 
termination.  Some  related  results  which  were  computed  previously  [9] 

are  also  included  in  Figure  49  for  completeness. 

Comparisons  of  calculated  and  measured  backscattering  results  for 

the  circular  inlet  case  are  considered  next.  The  measurements  were 
performed  at  the  Ohio  State  University  ElectroScience  Laboratory  [211. 
The  circular  inlet  may  be  terminated  inside  by  one  of  the  following 
structures:  a  planar  conductor  (short  circuit),  a  planar  impedance 
surface,  a  blade-disk  structure,  and  conical  and  hemispherical  hubs 
which  may  sit  on  a  planar  disk  or  a  disk-blade  structure.  In  Figures  50 
through  57,  the  solid  curves  are  the  measured  and  the  dashed  curves  the 
calculated  backscatteri ng  results.  Figures  50  and  51  show  the 
comparison  of  the  calculated  and  measured  backscattering  patterns  for  a 
circular  inlet  with  a  short  circuit  termination  at  three  different 

A  A 

frequencies  for  vertical  (<}>)  and  horizontal  (0)  polarizations, 
respectively.  Additional  horizontal  polarization  patterns  at  10  GHz  and 
10.2  GHz  are  shown  in  Figure  52  which  illustrates  the  sensitivity  in  the 
change  of  horizontal  polarization  patterns  with  frequency.  A  comparison 
of  backscattering  patterns  for  a  disk-blade  termination  is  shown  in 
Figure  53;  likewise,  they  are  shown  in  Figures  54  and  55  for  the  case  of 
conical  and  hemispherical  hubs  on  a  planar  disk,  respectively,  and  in 
Figures  56  and  57  for  conical  and  hemispherical  hubs  on  a  blade 
structure,  respectively.  It  is  observed  that  the  P.0,  approximation  for 
determining  the  reflection  from  the  termination  seems  to  give  a  good 
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overall  prediction  of  the  patterns,  except  in  a  few  cases  involving 
blade  terminations  where  the  patterns  may  not  compare  sufficiently  well 
with  measurements,  especially  for  the  horizontal  polarization  case  which 
appears  to  be  more  sensitive  to  parameter  changes  than  the  vertical 
polarization.  This  suggests  that  due  to  the  presence  of  the  blades,  the 
interactions  between  the  blades  themselves,  which  are  not  accounted  for 
in  the  P.0,  approximation,  could  possibly  have  a  non-negligible 
contribution  to  the  scattering  matrix  S .  Therefore,  a  refined  analysis 
could  be  used  to  improve  the  accuracy  of  the  calculated  patterns  for 
this  case.  In  addition,  the  experimental  model  for  the  blades  should  be 
more  carefully  made  than  at  present. 

Due  to  the  fast  varying  nature  of  the  RCS  patterns  with  frequency, 
some  frequency  scan  patterns  are  also  included.  In  particular,  results 
for  the  calculated  and  measured  frequency  scan  patterns  are  given  in 
Figures  58  through  61.  It  is  noted  that  the  calculated  backscattering 
returns  in  Figure  58  and  59  generally  are  stronger  than  the  measured 
ones.  This  is  due  to  the  fact  that  in  the  measurements,  the  terminal 
conducting  plane  does  not  completely  fill  up  the  inlet  cross  section, 
allowing  electromagnetic  energy  to  pass  through  and  this  results  in 
lower  measured  backscattering  returns.  It  is  important  to  note  that  the 
theoretical  predictions  show  the  same  general  trends  in  the  variation  of 
RCS  with  frequency  as  those  predicted  by  the  measurements. 
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Figure  51.  Comparisons  of  calculated  and  measured  horizontal 
polarization  backscatterlng  patterns  at  (a)  9  GHz, 
(b)  10  GHz,  and  (c)  11  GHz. 
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Figure  52.  Comparisons  of  calculated  and  measured  horizontal 
polarization  backscatterlng  patterns  at  (a)  10  GHz, 
(b)  10.2  GHz. 
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Figure  53.  Comparison  of  calculated  and  measured  backscatterlng 

patterns  at  10  GHz  (a)  vertical  polarization,  (b)  horizontal 
polarization. 
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Figure  53.  Continued. 
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Figure  54.  Comparison  of  calculated  and  measured  backscatterlng 
patterns  at  10  GHz,  (a)  vertical  polarization,  (b) 
horizontal  polarization. 
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Figure  55.  Comparison  of  calculated  and  measured  backscatterlng 
patterns  at  10  GHz.  (a)  vertical  polarization,  (b) 
horizontal  polarization. 
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55.  Continued. 
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Figure  56.  Comparisons  of  calculated  and  measured  backscatterlng 
patterns  at  10  GHz.  (a)  vertical  polarization,  (b) 
horizontal  polarization. 


114 


MAGNITUDE  IN  OB 

-HO.  -35.  -30.  -35.  -20.  -IS.  -10 


I 

I 

I 

I 


(b) 


Figure  56,  Continued, 
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Figure  57.  Comparisons  of  calculated  and  measured  backscatterlng 
patterns  at  10  GHz.  (a)  vertical  polarization,  (b) 
horizontal  polarization. 
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MAGNITUDE  IN  DB  MAGNITUDE  IN  DB 


(a) 


Figure  58.  Horizontal  polarization  backscatterlng  frequency  scans  with 
0l*O.  (a)  calculated,  (b)  measured. 
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MAGNITUDE  IN  DB  MAGNITUDE  IN  DB 


An  interesting  observation  as  a  consequence  of  using  the  P.0, 
approximation  in  determining  the  modal  reflection  coefficients  from  a 
conical  hub  termination  in  a  circular  inlet  is  depicted  in  Figure  62. 

As  shown  in  the  figure,  an  incident  mode  with  a  modal  ray  angle  B  is 
converted  by  the  conical  hub  termination  into  reflected  modes  which  are 
strong  if  their  modal  ray  angles  are  close  to  ir-2a-B  or  w-2a+B.  The 
ability  to  convert  an  incident  mode  into  several  reflected  modes  can 
have  a  drastic  effect  on  the  backscattering  return  when  the  terminal 
position  is  changed.  An  example  is  shown  in  Figure  63,  which  compares 
calculated  backscattering  returns  of  a  circular  inlet  with  a  conical  hub 
termination  at  different  positions.  A  simple  explanation  for  the  large 
change  in  the  backscattering  return  for  the  axial  incidence  (9=0)  case 
is  as  follows.  A  strong  TEn  mode  is  coupled  through  the  open  end  into 
the  circular  inlet.  Subsequently,  it  is  converted  by  the  concal  hub 
termination  into  reflected  TEn,  TE^  and  TE13  modes.  The  three  modes 
radiate  approximately  equally  in  the  axial  direction.  Thus  when 
changing  the  terminal  position,  although  the  relative  strength:  of  the 
three  modes  remain  virtually  unchanged,  their  different  phase 
combinations  can  result  in  significant  changes  in  the  backscattering 
returns  as  shown  in  Figure  63.  The  same  argument  can  also  be  applied  to 
the  non-axial  incidence  case,  although  the  situation  is  more  complicated 
there. 

If  the  base  length  (hd)  of  the  conical  hub  termiation  in  Figure  63 
is  reduced,  but  with  everything  else  remaining  the  same,  the 
backscattering  return  then  becomes  less  sensitive  to  the  change  in  the 


122 


terminal  position  of  the  conical  hub,  especially  in  the  axial  direction 
as  shown  in  Figure  64.  This  is  because  for  the  case  of  axial  incidence, 
only  the  TEn  mode  plays  an  important  role  in  contributing  to  the 
backscatteri ng  in  contrast  to  the  previous  case  of  Figure  63  wherein 
three  modes  contributed  equally  to  the  on-axis  backscatteri ng  return. 
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ir  -2.  a  +  f3 


Incidentand  reflected  modal  rays  associated  with 
hub  termination. 


a  conical 
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8  (DEGREES) 

(a)  vertical  polarization 

Figure  63.  Comparison  of  the  backscatteri ng  patterns  for  a  circular 
waveguide  with  conical  hub  terminations  at  different 
terminal  positions. 
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(b)  horizontal  polarization 


Figure  63.  Continued. 
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(a)  vertical  polarization 

Figure  64.  Comparison  of  the  backscattering  patterns  for  a  circular 

waveguide  with  conical  hub  termination  at  different  terminal 
positions. 
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FF.  PATTERN  (dB) 
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Figure  64.  Continued. 
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The  method  developed  in  the  previous  sections  can  also  be  applied 
to  evaluating  fields  scattered  by  open-ended  rectangular  and  semi¬ 
circular  inlets  mounted  on  a  ground  plane.  This  problem  can  be 
simplified  using  image  theory.  Thus,  this  problem  simply  reduces  to 
that  considered  before;  i.e.,  to  evaluate  the  fields  scattered  by 
open-ended  rectangular  and  circular  inlets  which  are  illuminated  by  two 
sources  correspondi ng  to  the  original  one  and  its  image. 

Figure  65  illustrates  a  comparison  of  the  calculated  backscattering 
results  from  an  open-ended  circular  waveguide  and  an  open-ended  semi¬ 
circular  one  mounted  on  a  ground  plane,  both  of  which  are  terminated  in 
the  interior  by  a  planar  short  circuit.  For  axial  incidence,  the 
existence  of  the  ground  plane  effectively  increases  the  backscattering 
return  by  6  dB  if  the  electric  field  is  polarized  perpendicular  to  the 
ground  plane,  and  it  reduces  the  backscattering  return  to  zero  if  the 
electric  field  is  polarized  parellel  to  the  ground  plane,  as  shown  in 
Figure  65.  Finally,  a  comparison  of  the  calculated  near  zone  bistatic 
scattering  from  an  open-ended  rectangular  waveguide  and  another  one  of 
half  the  size  mounted  on  a  ground  plane,  both  of  which  possess  a  short 
circuit  termination  inside,  is  shown  in  Figure  66.  It  is  interesting  to 

a 

note  that  for  the  <J>  polarization  case  (Figure  66b),  the  direction  of 
incidence  and  the  size  of  the  waveguide  are  chosen  in  such  a  way  that 
the  rim  scattering  of  the  incident  plane  wave  and  its  image  by  the 
rectangular  waveguide  mounted  on  the  ground  plane  exactly  cancel  each 
other.  Thus,  for  this  case,  the  entire  scattering  (dashed  curve)  is  due 
only  to  the  reradiation  of  the  field  reflected  back  by  the  short  circuit 
termination.  On  the  other  hand,  such  a  cancellation  of  the  rim 
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FF.  PATTERN  (dB) 


(a)  0  polarization 

Figure  65.  Comparison  of  backscattering  returns  from  an  open-ended 
circular  waveguide  (solid  curve)  and  an  open-ended 
semi-circular  waveguide  mounted  on  a  ground  plane  (broken 
curve). 
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(b)  <j>  polarization 


Figure  65.  Continued. 
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(a)  0  polarization 

Figure  66.  Comparison  of  near  zone  bistatic  scattering  from  an 

open-ended  rectangular  waveguide  (solid  curve)  and  another 
one  of  half  the  vertical  dimension  mounted  on  a  ground  plane 
(broken  curve)  for  observation  at  (r*4X,<|>*90°). 
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scattering  by  the  image  may  not  occur  for  the  0  polarization,  because 
the  main  contribution  in  that  case  comes  from  the  vertical  plates. 

V.  CONCLUSION 

The  work  described  in  this  report  provides  relatively  efficient 
solutions  which  are  based  on  a  combination  of  high  frequency  and  modal 
techniques  for  the  problems  of  near  and  tar  zone  bi static  scattering  of 
EM  plane  waves  incident  on  open  ended  semi -infinite,  perfectly 
conducting  circular  and  rectangular  inlet  configurations  with  an 
interior  termination.  The  interior  terminations  could  be  either  a 
planar  perfectly  conducting  surface  *r  a  planar  impedance  (or 
dielectric)  surface  termination.  It  is  noted  that  when  the  impedance 
vanishes  then  the  impedance  surface  termination  reduces  to  the 
perfectly-conducting  case. 

Some  additional  simple  interior  terminations  have  also  been 
considered  in  this  work  for  the  case  of  the  circular  inlet  geometry; 
they  are: 

a)  a  conducting  disk -blade  termination 

b)  a  conducting  conical  or  hemi -spherical  hub  placed  on  a 
conducting  disk,  or  on  a  conducting  disk-blade  structure. 

Moreover,  the  related  problems  of  the  EM  scattering  by  an  open 
ended,  semi -inf inite  perfectly  conducting  semi-circular  and  rectangular 
inlets  on  a  perfectly  conducting  ground  plane  have  been  analyzed 
directly  via  an  application  of  image  theory  In  combination  with  the 
solutions  which  have  been  obtained  for  the  circular  and  rectangular 
Inlets  without  the  ground  plane. 


It  is  seen  that  the  RCS  of  the  above  structures  are  quite  sensitive 
to  frequency.  Also,  the  scattering  patterns  of  the  horizontally 
polarized  wave  are  very  sensitive  to  small  changes  in  structure.  The 
calculations  for  the  far  zone  EM  backscatteri ng  case  which  have  been 
performed  in  this  work  are  compared  with  corresponding  measurements. 

The  overall  agreement  between  the  calculation  and  measurement  is  quite 
good.  No  measurements  have  been  performed  on  this  contract  for  the  case 
of  the  near  zone  bistatic  scattering  by  such  inlets.  In  the  future,  it 
would  be  worthwhile  to  make  near  and  far  zone  measurements  of  scattering 
by  somewhat  more  improved  experimental  models  than  those  used  presently. 
Also,  it  would  be  worth  extending  the  research  effort  to  find  ways  to 
analyze  the  more  complex  problem  of  scattering  by  non-uniform  inlet 
geometries  in  a  relatively  efficient  manner. 
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APPENDIX  I 


WAVEGUIDE  EXCITATION  PROBLEM 


Let  a  source  S  be  located  inside  waveguide  as  shown  in  Figure  A.l. 
S  can  be  an  electric  current  source  J  or  a  magnetic  current  source  M, 
which  generates  E*,  H*.  Thus, 


s±  ‘l  ap  (Ep-  "P 


(A.l) 


—  A 

are  the  fields  generated  by  S  in  the  ±z  directions.  In  this  summation 
p  is  the  compact  summation  index  representing  the  double  summation  over 
mode  indices  "n"  and  "m".  It  is  of  interest  to  find  (E*;  fC).  Here, 

E±  will  be  determined  via  an  application  of  the  reciprocity  theorem  as 
In  [13]. 

The  modal  fields  can  be  decomposed  into  transverse  and  axial 
components  as  follows: 


(A. 2a) 

Ep  *  ApC£pt  * 

(A. 2b) 
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H.  =  h  ;  H  =0  ;  E  =  i  and  E  =  e 

Pt  nm  pz  pt  nm  pz  z,nm 

Case  (1): 

It  is  of  interest  to  find  the  strength  of  the  mode  E",  H" 

q  q 

generated  by  S.  (Here,  q  is  the  compact  index  representing  the  mode 
indices  "ij".) 
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From  the  reciprocity  theorem: 


///  E+  •  J  dv  for  S  =  J 
V  q 

-///  H+  •  M  dv  for  S  =  M 
V  q 

(A. 3) 

where  E,  H  =  E+,  H+  on  S+  and  E,  H  =  E  ,  H  onS.  Also,  E  *  n  =0 
±  -  q  Sw 

and  E0  x  n  =  °* 

SW 

Thus,  the  above  reciprocity  relation  simplifies  to: 


j  j  [E  x  H+  -  E+  x  h]  •  n  ds  = 

S  s  +S  q  q 

V-  W 


//  [E+  x  H+  -  E+  x  H+]  •  2  ds  +  JJ  [E“  X  ii+  -  E+  x  H"]  *(-z)ds 

S  M  H  $  H  H 


III  E*  •  J  dv 

V  q 

-III  H*  •  M  dv 

V  q 


for  S  =  J 

for  S  =  M 


(A. 4) 


Substituting  E+,  H+  and  E",  H“  from  Equation  (A.l)  and  using  the 
orthogonality  condition: 

JJ  x  H*]  •  n  ds  =  0  for  q  *  p  (A. 5) 

C  r  H 

one  obtains  a  relationship  in  terms  of  only  qth  mode  (i.e.,  p=q  case) 
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ff  a+  ffE  +  E  )  x  (H  .  +  H  )e 
q  u  qt  qzJ  1  qt  qzJ 


“J2V+ 


-j2B  z 

(Eqt  *  E,z>  *  (Hqt  +  Hqz>e  1  ’  2  dS 


ff  A"  ffE  -  E  )  x  (H  .  +  H  ) 
JJ  q  u  qt  qzJ  ^  qt  qz; 


-fE  +  E  )*(-H+H  )]  •  z  ds 
lqt  qzJ  ^  qt  qz;j 


///  E+  •  5  dv  for  S  =  5 
V  q 

-///  H+  •  M  dv  for  S  =  M 
V  q 


(A. 6) 


or; 


-2A"  ff  E  x  H  ■  z  ds 
q  qt  qt 


-  -jB  z 

///  (E  +  E  )  •  J  e  q  dv  for  S=J 
V  qt  qz 

-  -JB  z 

-///  (H  +  H  )  •  M  e  q  dv  for  S=M 
y  qt  qz 

(A. 7) 


Thus, 


A"  = 

q 


2<S  Eqt  *  Hqt 


z  ds 


111  <Eqt  *  V  •  5  e  JV  d'  f°r  54 

-  “JB  z 

-fff  (H  *  +  H  )  •  He  q  dv  for  S=M 
v  qt  qz' 

(A. 8) 
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where  A~  is  the  excitation  coefficient  of  the  mode  E"  traveling  to  the 
q  q  3 

right  of  the  source,  i.e.,  in  the  region  z<z-  shown  in  Figure  A.l. 


Case  (ii): 

If  E”,  H"  is  chosen  as  the  waveguide  mode  of  interest  as  in 

Figure  A. 2,  then  the  coefficient  A+  of  the  mode  E+  traveling  to  the 

q  q 

left  of  the  source  (z>z+)  is  given  by: 


A+  = 


2/J  E  x  H  . 
s  qt  qt 


z  ds 


it,i  -  V  • 3  *°V 

-  j3„z 

-///  (-H  .  +  H  )  -He  q  dv 
V  qt  qz 


(A. 9) 


Note  that  //  Eqt  *  «qt^  ds  -  //^  x  ds  -  //  Eqt  x  ds 

where  Sc  is  any  cross-sectional  area  of  the  waveguide. 

Although  Equations  (A. 8)  and  (A. 9)  are  given  for  a  volume  source 
distribution,  it  is  very  simple  to  modify  them  for  a  surface  (or  line) 
source  distribution.  For  this  purpose,  the  volume  integral  on  the 
right  hand  side  of  Equations  (A. 8)  or  (A. 9)  should  be  replaced  by  a 
surface  (or  line)  integral  over  the  extent  of  the  surface  (or  line) 
source  distribution. 
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Figure  A. 2.  Waveguide  geometry  with  fields  radiated  to  the  left  (z>z+) 
of  being  of  Interest. 
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APPENDIX  II 


TABLES  FOR  THE  CIRCULAR  AND  RECTANGULAR 
WAVEGUIDE  MODE  FUNCTIONS 
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TABLE  1 


PROPERTIES  OF  MODES  IN  CIRCULAR  WAVEGUIDE 
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TABLE  2 

PROPERTIES  OF  MODES  IN  RECTANGULAR  WAVEGUIDE 


TE. 


TM 


2  2 
na  +mb 

h  =  N  Y  “7T - cos  n  x  cos  ray 

z  nm  nm  J0nm  a  d 


ez  *  0 
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nm 
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V 
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APPENDIX  III 


EXPLICIT  RADIATION  COEFFICIENT 
FOR  CIRCULAR  AND  RECTAN6ULAR  INLETS 

l"W 

As  mentioned  In  the  text,  the  coupling  coefficient  [S21]  is 
related  to  the  radiation  coefficient  [S12]  by  reciprocity. 

Therefore,  if  one  is  known,  the  other  can  be  easily  calculated 
from  reciprocity  relationship  given  in  Section  III.D. 

For  the  near  zone  radiation  case,  there  is  no  simple  closed  form 
analytic  solution  that  can  be  obtained  to  calculate  these  coefficients; 
whereas,  in  the  far  zone  radiation  case,  one  can  find  simple  closed 
analytic  forms  for  the  radiation  coefficients  and/or  coupling 
coefficients. 

In  the  following,  the  explicit  form  of  the  radiation  coefficients 
in  the  far  zone  are  given  for  both  circular  and  rectangular  waveguide 
modes . 

A.  For  Circular  Inlets: 

In  this  case,  the  far  zone  field  can  be  written  as 


where  Eq  and  E^  can  be  separated  as  contributions  from  the  Klrchhoff 
approximation  and  the  equivalent  Ufimtsev  edge  currents. 
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Ee  =  Eek  +  E0u 


ea  =  ea.  +  ea 

$  <j>k  <j>u 


The  contribution  from  the  Kirchhoff  approximation  becomes: 


TERm  modes  incidence: 

1+COS0  C0S6 


Eq.  =  jnkZ  N  n 
0k  d  o  nm 


nm 


W  Jn  (kas(n9) 


sin6 


nm 


E ..  =  jnkZ  N  P  - - - 7T  J  (P  )  J  (kasina) 

<j>k  o  nm  nm  2(cos5nm-cos0)  n  nm  n 


6  =  cos'1 (3  /k) 

nm  nm 


modes  incidence: 


E0k  ”  J  kNnm  Pnm  2(cos5nm-cos0)  unvrnm/  un 

E..  =  0 
<t>k 


sin0 


J  (P  J  J  (kasinO) 


Likewise,  the  contribution  from  the  equivalent  Ufimtsev  edge  currents  is 
given  by: 

modes  Incidence: 

0 

Efi  =  jnZA  N  nf (0,5  )  J  (p‘  )  [enm  slnT  j"  (kas1n0) 

0u  o  nm  nm  n  nm  nm  c.  n 

6. 


nm  cos0 
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The  contribution  from  the  Klrchhoff  approximation  becomes: 
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Likewise  the  contribution  from  the  equivalent  Ufimtsev  Edge  Currents  is 
given  by: 
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where 


u,  v,  A+,  B+  are  defined  as  before,  and 
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APPENDIX  IV 


EXPLICIT  MODAL  REFLECTION  COEFFICIENTS 
AT  THE  OPEN  END 

FOR  THE  CIRCULAR  AND  RECTANGULAR  INLETS 

As  mentioned  in  the  text,  the  modal  reflection  coefficients  at  the 
open  end  can  be  calculated  from  the  equivalent  magnetic  currents  and 
M*d,  i.e., 


/(hpg-M*  +  h2pq«M)da 

%q;nm  =  2;/?  x  *  *  ds 

pq  pq 

and 
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So  that 


pq;nm 


1  /  (h  -M* 
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+  h  •M*’d)dJl 
zpq 


The  advantage  of  using  this  form  Is  that  Rpq.nm  satisfies  the 
reciprocity  principle,  i.e., 


pq;nm  nm;pq 
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An  explicit  expression  for  the  reflection  coefficient  from  an  nm 
incident  mode  to  pq  reflected  mode  can  be  found  from  the  explicit  form 
of  above  equations,  which  is  given  by: 


A.  For  Circular  Inlets: 
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-j(2ka-ir/4) 

A  =  (-l)n  - 

4/ iika 


c  =  1  +  I  (jB)V  /i+r 
£=1 


D  =  1  +  l  ( —  jB )  /  a TT 
£  =  1 


B  -  (-!)"?  e'J'2ka 


Due  to  the  orthogonality  of  the  trigonometric  functions,  cosn<|>  and 
sinn<t>,  which  appear  in  the  model  functions,  the  above  reflection 
coefficients  apply  only  to  modes  with  the  same  <f>  variation  and  there  is 
no  coupling  between  modes  with  different  4>  variations.  It  is  noted  that 
the  effect  of  all  the  multiply  interacting  rays  across  the  aperture  are 
included  in  the  above  expressions  for  the  circular  inlet  opening;  in 
particular,  that  information  is  contained  in  the  "f"  and  "g"  terms  given 
above.  Unfortunately  it  Is  not  a  simple  matter  to  include  multiple 
interactions  for  the  rectangular  inlet  opening  which  is  dealt  with 
next. 
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B. 


For  Rectangular  Inlets: 
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N  ,  Y  „ .  n  .  and  m.  are  defined  as  In  Table  2,  Appendix  III, 
nm  nm  a  d 
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APPENDIX  V 


DEFINITIONS  OF  THE  GENERALIZED  SCATTERING  MATRICES 
FOR  THE  FAR-ZONE  SCATTERING  SITUATION 


When  the  observation  point  receeds  to  Infinity,  i.e.,  the  far  zone 
case,  the  scattering  matrices  are  defined  as  follows  [9]: 

[Sn]  Matrix: 


[S12]  Matrix: 
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And  the  definitions  of  [S21 3 >  [S22]  and  [Sr]  remain  the  same  as 
before.  Note  that  in  this  far  zone  case  the  range  (r)  dependence  is 
brought  out  of  [S^]  and  [S12L  and  the  entry  of  the  matrices  has  been 
changed  from  x,y,z  components  to  6,<fr  components. 


APPENDIX  VI 


NODAL  REFLECTION  COEFFICIENTS  FOR  DISK-BLADE 
AND  HUB  STRUCTURE  TERMINATIONS  IN  CIRCULAR  INLETS 


To  find  the  reflection  coefficients  from  the  disk,  disk-blade  and 
hub  structures  in  circular  waveguides  the  method  of  P.0,  approximation 
is  employed.  Although  the  P.0,  approximation  does  not  take  into  account 
the  proper  boundary  conditions  at  the  end  of  the  structure,  it  certainly 
picks  up  the  dominant  effect  and  yields  good  approximation  to  the 
scattering  results  as  compared  with  the  measurements. 

Due  to  the  symmetry  in  the  geometries  of  the  terminal  structures 
studied  here,  the  coupling  coefficients  between  modes  with  different  <|> 
variations  are  small,  so  that  they  are  neglected  in  the  following 
discussion.  Only  the  mutual  coupling  between  those  modes  with  the  same 
<J>  variation  are  Included. 

An  explicit  reflection  coefflcientt  for  a  disk-blade  termination 
is  given  by 
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where  the  Integration  from  p=0  to  p=b  is  over  the  disk  region  and  that 
from  p*b  to  p»b+d  is  over  the  blade  structure;  a  is  the  angle  extension 
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of  the  blade  and  0  the  angle  extension  of  the  air  slot;  Z  =  Zq  for 

a  TE  mode  and  Z  for  a  TM  mode.  The  model  functions  are  those 
nm  k  o  nm 

given  in  Appendix  II. 

For  a  hemispherical  hub  termination  with  radius  b  the  reflection 
coefficients  are 
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Similarly,  for  a  conical  hub  termination  with  base  2b  and  height  d,  the 
reflection  coefficients  are 
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